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Introduction 



Obtaining information on short range correlations (SRC) in nuclei is a pri- 
mary goal of modern nuclear physics pQ. Interest in SRC stems not only from 
the necessity to firmly establish the limits of validity of the standard model of 
nuclei, i.e. a non relativistic description in terms of two- and three-nucleon 
interactions, but also from the impact that the knowledge of the detailed 
mechanism of SRC would have in understanding the role played by quark 
degrees of freedom in hadronic matter and the properties of the latter in 
dense configurations [2] . Recently, evidence of SRC has been provided by new 
experimental data on inclusive [A{e, e')X] [3l 0] and exclusive [A{e, e'pN)X 
and A{p,pN)X] lepton and hadron scattering off nuclei at high momentum 
transfer (Q^ > 1 GeV^) (see Ref. [3] and references therein quoted). In inclu- 
sive scattering the observation of a scaling behavior of the ratio of the cross 
section on heavy nuclei to that on the Deuteron [3], for values of the Bjorken 
scaling variable 1.4 < xbj ^ 2, and to that on ^He |4j, for 2 < < 3, has 
been interpreted as evidence that the electron probes two- and three-nucleon 
correlations in complex nuclei similar to the ones occurring in the two- and 
three-nucleon systems P, [7] . By combining the results of inclusive and exclu- 
sive experiments, a convincing experimental evidence of SRC in nuclei have 
been eventually found. It should be pointed out, however, that in inclusive 
experiments statistics in the region of three-nucleon short range correlations 
is very poor and, at the same time, a general framework to describe these 
correlations is still lacking; moreover, in exclusive experiments, the found ev- 
idence of SRC is limited to the ^^C nucleus. For these reasons it is not only 
necessary to extend experimental measurements to other nuclei, as planned 
e.g. at the Thomas Jefferson National Accelerator Facility (JLab), but it is 
also urgent to improve our theoretical knowledge on the nature of two- and 
three-nucleon correlations. The aim of this Thesis is to critically review the 
present theoretical and experimental knowledges on SRC in nuclei, and, at 
the same time, to provide a theoretical framework within which to coherently 
treat two- and three-nucleon correlations. As already pointed out, most of 
our experimental knowledge on two- and three-nucleon correlations comes 
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from exclusive and inclusive experiments on lepton and hadron scattering 
off nuclei at high momentum transfer. It should be pointed out, however, 
that whereas exclusive processes can directly access the relative and center of 
mass motions of a correlated pair in a nucleus [HlIHllIO], obtaining information 
on these quantities from inclusive scattering is, in principle, more difficult. 
Various approaches based on scaling concepts have therefore been proposed, 
going from the scaling behavior of the cross section ratio plotted versus xbj, 
to the scaling behavior of the ratio of the nuclear to the nucleon cross sections 
plotted versus proper scaling variables; among the latter, a process that has 
been most investigated in the past is the so called Y-scaling, for it is believed 
that this may represent a powerful tool to extract the high momentum part 
of the nucleon momentum distribution which is governed by SRC fHX \T2l . 
In this Thesis the concepts of F-scaling will be critically reanalyzed, mainly 
because of: i) the lack of a general consensus about the usefulness of such 
a concept, and ii) a strong renewal of interest in Y-scaling owing to recent 
experimental data on A{e,e')X reactions from the JLab O |T2]. We will 
show that the analysis of inclusive scattering in terms of proper Y-scaling 
variables could indeed provide useful information on SRC; to this end, fol- 
lowing the suggestion of Refs. [HI [TTl HHj, a new approach to Y-scaling and 
its usefulness will be illustrated in detail. 

In the first part of the Thesis, mainly in Chapter 1, we present the necessary 
formalism to introduce the original part, which is illustrated in Chapter 2-6. 
In more detail, the structure is as follows: 

1. In Chapter 1, the realistic many-body problem of nuclei and cold 
hadronic matter is recalled, in order to introduce the concept of SRC; 
an overview of the experimental evidence of SRC is given, and the 
relevance of SRC in various fields of Physics is illustrated. 

2. In Chapter 2, the basic features of the spectral function and the nucleon 
momentum distributions, which are two basic quantities of our new 
approach to inclusive cross sections, are recalled. 

3. In Chapter 3, the formalism of inclusive electron scattering off nuclei 
and the general expression of the cross section is presented, within the 
plane wave impulse approximation (PWIA), and by taking into account 
the final state interaction (FSI) of the knocked out nucleon with the 
residual system {A — 1). 

4. In Chapter 4, the inclusive process is analyzed in terms of y-scaling: 
three different scaling variables are introduced, each one describing a 
particular process occurring in electron scattering off nuclei. 



V 



5. In Chapter 5, the results of our calculations of the xbj dependence 
of the inclusive cross section ratios in terms of PWIA and FSI are 
presented. 
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Chapter 1 

Short range correlations in 
nuclei 

Introduction 

We will now introduce what short range correlations are and what is their 
relevance in modern physics. 

1.1 The realistic many-body problem of nu- 
clei and cold hadronic matter 

In what follows, we will consider a bound system of Z protons and neu- 
trons, with A = Z ^ N . We will simply call such a system a nucleus A. 
What exposed in the following equally well applies to both nuclei and cold 
hadronic matter, the latter being a system composed of an infinite number 
of bound nucleons (e.g. neutron stars). As stressed in Ref. [19], when par- 
ticles interact with each other through the intervening mechanism of a field, 
the description of their dynamical behavior in terms of instantaneous poten- 
tials is only an approximate nature, and the two-body V2{xiiXj\ three-body 
t'3(xi, Xj, Xfc),. . . , A-body potentials may be regarded as successive stages of 
this approximation. Thus, considering the nucleus A as a non relativistic 
quantum-mechanical system, its quantum states are the solution of the nu- 
clear many-body problem, represented by the following Schrddinger equation 

H[xx . . . xaWa[xx ...xa) = En^\{xx ...xa) (1.1) 
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with Hamiltonian 

j_2 A A A 

H{Xi . . .Xa) = y^^t^+y^ V2{Xi,Xj)+ V3{Xi,Xj,Xk) + . . . 

1=1 i<j=l i<j<k=l 

(1.2) 

Here Xi = {rj, Sj, tj} denotes the generalized coordinate of the i — th nucleon, 
which includes its radial coordinate rj, spin Sj and isospin tj; n stands for the 
set of quantum numbers of the state under consideration; m^r is the nucleon 
mass. 

Solving the nuclear many-body problem is not an easy task; the reasons are 
many fold: 

• many-body forces are unknown; 

• the two-body potential obtained from the analysis of nucleon-nucleon 
(NN) scattering data is very complicated, owing to its spin, isospin and 
tensor dependences [20] : 

• it is not yet clear what is the role (if any) played by the quark-gluon 
structure of the nucleon in the description of nuclear properties; 

• it is not yet clear to which extent the nucleons bound in a nucleus 
retain the same properties as the free ones; 

• last but not least: can nucleon motion in a nucleus be considered within 
the non relativistic approximation? 

It has however been demonstrated in Ref. [19] that, independently of the 
detail of the field, the m-body potentials, in systems governed by the strong 
force, can be written as follows 

(m-body potentials) ~ y — j x (two-body potentials) (1.3) 

where m = {3, . . . , A}, and is the average nucleon velocity which, using 
the uncertainty principle and the known nuclear dimension, can be estimated 
to be vn ~ 0.02c, which means that, to a large extent, nuclear systems 
can be considered as non relativistic systems bound by two- and, at most, 
three-nucleon interactions. Such a conclusion is confirmed by a wealth of 
experimental information on basic properties of nuclei (e.g. the dependence of 
their radii upon A^/^, which, in turns, leads to a constant value of the binding 
energy per nucleon and to the Constance of the volume, the similarity between 
the magnetic moments of odd nuclei and the nucleon magnetic moments, etc.) 
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hadronic matter 



which lead to the conclusion that the atomic nucleus can, on the average, 
be described as an incompressible low density system (with density po ~ 
0.17 N/ fm^) composed of non relativistic nucleons interacting mainly via the 
same two-nucleon strong force acting between free nucleons plus, at most, we 
reiterate, three-nucleon forces. The nuclear many-body problem thus reduces 
to what has been called the standard model of nuclei [21] , described by the 
following Schrddinger equation 



2m 



N 



A A 

i<j=l 



1=1 



i<j<k=l 



Kn(xi...x^)(1.4) 



Even in this simplified form, Eq. fll.4p is difficult to solve, due to the com- 
plicated structure of the two-nucleon interaction. For such a reason, in the 
past half century, various phenomenological models have been proposed to 
explain the structure of nuclei, and the nuclear shell model (SM), for which 
Maria Goppert-Mayer and Hans D. Jensen were awarded by the Nobel Prize 
in 1963, turned out to be the most successful one [22]. In the simplest version 
of this model, the independent particle shell model, the nucleus is described 
as an ensemble of independent nucleons which move in an average potential 
filling, according to the Pauli Exclusion Principle, proper shell model states. 
Moreover, neutrons and protons occupy all states below the Fermi level, leav- 
ing the above states empty. More technically, one says that the occupation 
probability of states below the Fermi level is one, and above the Fermi level 
is zero. 

The Hamiltonian of the system reduces to 



Ho{xi ...xa) 



and the Schrodinger equation is 



2m 



N 



;i.5) 



2m 



N 



5o = eo <Po 



1.6) 



$0 being a Slater determinant. In the most refined SM description, various 
types of residual interactions are added to Eq. (11. 5p : these include, for ex- 
ample, the spin orbit interaction, non spherical single particle potentials to 
account for the deviation from the spherical shape of classes of nuclei in the 
periodic table, and others. The success of the advanced SM in reproducing 
many properties of nuclei was awarded in 1975 by a Nobel Prize to A. Bohr, 
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B. Mottelson and L. Rainwater. In the advanced shell model, a prominent 
role is played by the so-called long range correlations, whose main effect is 
to partly deplete the occupation probability of the states below the Fermi 
level, making the states above the Fermi level partially occupied. It should 
however be pointed out that the main feature of the independent particle 
and advanced shell models, is the independent particle motion. This fact is 
difficult to reconcile with one of the main features of the realistic NN inter- 
action, namely the strong repulsive core at relative distances of the order of 
0.5 — 0.6 fm, which is one of the facts which makes the solution of the nuclear 
many-body problem (11. 4p a very difficult one. 

The realistic two-body interaction, which explains two-body bound and scat- 
tering data, has the following form [20] 

N 

v,{x,,x^ = Y,v^ir,^0^ (1-7) 

n=l 

where r^j = |rj — r^l is the relative distance of nucleons i and j, and n, 
ranging up to = 18, labels the state-dependent operator 0^"^ whose first 
six components are defined as follows 





= 6^. = 1 


(1.8) 






(1.9) 






(1.10) 






(1.11) 




= O^*^ - 9-- 


(1.12) 


1-3 




(1.13) 



Here 

Sij = 3{(Ti- fij) {(Tj ■ fij) - {(Ti ■ (Tj) (1.14) 

is the tensor operator, and cTj and Tj are the spin and isospin of the z-th 
nucleon of the pair, respectively. It can be seen that the two-body interaction 
exhibits a strong state dependence. Particular worth being mentioned are the 
strong short range repulsion and the tensor attractive interaction in T = 
and iS* = 1 states, as shown in Fig. II. 1[ These features of the NN interaction 
generate, as we shall see in more detail in the following, strong correlations 
between nucleons, which are not present in a shell model picture. 
Recently, important progress has been made in solving the many-body Schro- 
dinger equation (11.41) . which is rewritten by considering only the ground state. 
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Figure 1.1: The Paris nucleon-nucleon potential in triplet {S = 1) central 
and tensor states. After Ref. [23]. 
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denoting = and Eq = Ea, namely 



2m 



N 



A A A 

i=l i<j=l i<j<k=l 



= Ea'^a- 
(1.15) 

The most advanced approaches to the solution of the nuclear many-body 
problem rely on the numerical integration of Eq. fll.lSp by Monte Carlo 
techniques |23], in particular the Variational Monte Carlo (VMC) method, 
which is used to optimize the expectation value of observables by adjusting a 
trial wave function, and the Green Function Monte Carlo method. By Monte 
Carlo techniques it was possible to solve Eq. fll.lSp for both the ground and 
excited states in the range 3 < A < 8, as shown in Fig. 11.21 For heavier 
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Figure 1.2: Variational Monte Carlo (VMC) and Green Function Monte Carlo 
(GFMC) energies using the two-body AV18 [25] and three-body UIX [26] 
interactions compared with experiment (Exp). Black dashed lines show the 
indicated breakup thresholds for each method. The Monte Carlo statistical 
errors are shown by the light blue and yellow bands. After Ref. 



nuclei the Monte Carlo methods become very difficult to apply, and alterna- 
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tive realistic approaches have been developed, such as the Hypernetted chain 
method with correlated wave function [27], large shell model basis approach 
[28], the coupled-cluster theory [29], and the variational calculations with 
correlated wave functions [30]. In this Thesis we will consider the number 
conserving cluster expansion approach with correlated wave functions [21] • 
In this method, the variational principle is applied in minimizing the expec- 
tation value of the Hamiltonian 



<H>-- 



> En 



with the trial nuclear wave function cast in the following form 
where 0o is the shell model mean field wave function, and 

n ^"(^-^^■) 

i<j 

is the correlation operator; the latter is defined in terms of the symmetriza- 
tion operator S, and the two-body correlation function 



(1.16) 

(1.17) 
(1.18) 



N 



r-) 



(n) 



(1.19) 



n=l 



with the same operatorial dependence appearing in the two-body potential 
V2{xi,Xj) given by Eq. (11.71) . In Ref. [31], a new effective method for the 
calculation of the expectation value of any quantum-mechanical operator A 
in the many-body ground state described by the wave function \E'^, i.e. 



<A>-- 



(1.20) 



with given by Eq. f ll.l7p . based upon the cluster expansion, is presented. 
Introducing the quantity 

Vi, = fl k - 1 (1-21) 
the expectation value (11.201) can be written as follows 



<A> 



<0o|FtiF|0o> 
< 0o|0o > 



< 00 


ni<,(i+%) ^ 


00 > 


< 00 




00 > 



< 00 


+ Si<j Vij + Yl,(i<j)<{k<l) VijVkl + ■ ■ 




00 > 


1+ < 00 


Si<j Vij 


00 > + . . . 



;i.22) 
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Expanding in series the denominator 

1 



l + x 



1 — X + X 



:i.23) 



one gets 



<A> 



< 0o|i|0O > + < (pQ\^flijA\ 



> + 



X 



1- < 



i<j 



> + < 



i<j 



i<j 



(1.24) 



and collecting all terms containing the same number of function fjij, one 
obtains the infinite series 



< A >=< A>o + <A>i + <A>2+...+ <A>n + 
At second order in t], one has, explicitly, 

< i >o = < (j)o\A\(Po > 



< A>i 

<A>2 



< 



< 



^01 

i<j 

{i<j)<{k<l) 

X <0o|^%|0O> 

i<j 



> 



7]ij A\(j)Q > - < i >o < 0o| X^^ijl 

7]ij7]uA\(t)0 > - < 0o| y~''7ij^|0O > 



:i.25) 

;i.26) 
1.27) 



X 



< 0o| fiijfjki\(j)o > - < 4>o\^Vij\4>o >^ |(l-28) 

(i<j)<{k<l) i<j 



where the term of the order of n contains fjij (fij) up to the rath (2dn) 
power. In Ref. [31], the ground state properties of closed shell nuclei, '^He, 
^^O and ^^Ca, have been calculated by minimizing the ground state energy 
at 2nd order in fjij with respect to the single particle wave functions and the 
correlation functions f^^\rij). The results for the latter, corresponding to 
the V8' interaction |32], are shown in Fig. 11.31 versus the distance r, in case 
of ^^O. It should be pointed out that: 

• at relative distances significantly smaller than the average internucleon 
distance (~ 1.7 /m), nucleons feel the strong central repulsion and the 
tensor attraction; 
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1 2 3 4r[fm]5 



Figure 1.3: The correlation functions f ll.l9p for ^^O vs. the relative distance 
r = rij, corresponding to the full Argonne V8' interaction used in Ref. [3T] . 



• at large distances, the effects due to noncentral correlations vanish, the 
central component is equal to 1, and the mean field wave function in 
Eq. fll.l7p is recovered. 

The net effect generated by the correlation function is demonstrated in Fig. 
II. 4[ where the calculated two-body density distribution for p-n and p-p pairs 
in ^^O resulting from the solution of Eq. f ll.lSp with interaction (II. 7p . con- 
taining the Argonne f 14 two-nucleon pi] and Urbana VII three-nucleon po- 
tential j35j, is compared with the shell model density solution of Eq. (11. 6p . It 
can be seen that, as a result of the contrasting short range central repulsion 
and the intermediate range tensor attraction, the realistic two-body density 
strongly differs from the shell model density in the region < r < 1.1 fm; 
since the average internucleon distance is ~ 1.7 — 2 fm, we will call the 
deviation from the shell model in the range < r < 1.1 fm short range cor- 
relations (SRC), always remembering that they are due to both short range 
central repulsion and intermediate range tensor attraction, with the latter 
acting only in T = and S = 1 states, and thus enhancing the n-p distri- 
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Figure 1.4: The two-body density distribution versus the relative distance 
r = Vij. Dot-dashed line: p-p mean field (MF) contribution, calculated from 
0o; solid line: p-p correlated contribution, calculated from \E'^ = \E't,; dotted 
line: p-n MF part; dashed line: p-n correlated component. After Ref. [53] 



butions. In what follows, it will be shown that SRC strongly depopulate the 
states below the Fermi sea, creating highly excited virtual two-particle (2p)- 
two-hole {2h) states. 



1.2 Experimental investigation of long and 
short range correlations in nuclei 

A simple cartoon of NN SRC, as predicted by theoretical calculations, is 
depicted in Fig. II. 5[ and the question arises as to whether such a picture can 
be experimentally observed. The answer is a positive one, and will be given 
qualitatively in what follows. 
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Figure 1.5: A simple cartoon of NN short range correlations. When the 
distance between two nucleons becomes smaller than the average internucleon 
distance, of order ~ 1.7 /m, nucleons result to be in a correlated pair, and 
their local density becomes comparable with those in the core of neutron 
stars, up to ~ 7 times the average nuclear density po- 



1.2.1 Exclusive lepton scattering: the A{e,e'p)X reac- 
tion 

A way to learn about correlations in nuclei is represented by the nuclear 
reaction in which an energetic electron e, with energy ei, knocks out from 
a nucleus A a proton p, which is detected in coincidence with the scattered 
electron e', with energy €2- In this process, the [A — 1) nucleus is left in 
some excited states with energy E* = u — Tp — Ta-i, where = ei — 62 is 
the energy transfer (the energy lost by the electron in the scattering process) 
and Tp and T^-i are the kinetic energies of the proton and the residual 
nucleus, respectively. This kind of nuclear ionization experiments, initiated 
more than 40 years ago [36], depends upon the number of protons in various 
shell model states (see §2.3p or, in other words, on the occupation probability 
in various shell model states, which, as we know, in the single particle shell 
model is one, for states below the Fermi sea (e.g. lsi/2 and P3/2 in ^^C), and 
zero for states above the Fermi sea. 

The situation for the ^^C nucleus is illustrated in Fig. 11.61 It can be seen 
that the first experiments, performed at low energy resolution, identified only 
two shell model states in ^^C, namely Ip and Is, whereas the high energy 
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•-• 

•-• 1s 

SRC (2p2h) 

(c) 

Figure 1.6: Nuclear ionization experiments ^^C(e, e'p)X performed at: (a) 
Laboratori Nazionali di Frascati (LNF) at low energy resolution and low re- 
moval energy Emin < E < AOMeV [36j, (b) the National Institute for Nuclear 
Physics and High Energy Physics (NIKHEF) at high energy resolution and 
low removal energy -Emm < E < 2OMeV'07l and (c) the Jefferson Laboratory 
(JLab) facilities at high values of the removal energy Emin < E < 200 MeV] 
the black area is attributed to SRC 1381. 
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resolution experiments [3Z] demonstrated the occupation [30] of bothp3/2 and 
Pi/2 shells, as a result of long range correlations {configuration mixing [39]). 
The occupation numbers of the valence protons in various nuclei obtained 
from (e, e'p) reactions are summarized in Fig. 11.71 It can be seen that 
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Figure 1.7: Spectroscopic strength for knocked out valence protons measured 
with the reaction (e, e'p) relative to independent particle shell model predic- 
tion HO]. 

the occupation probability of the least bound ('valence') nucleons is about 
0.6 — 0.7; it strongly deviates from one, i.e. from the prediction of the 
independent particle shell model. Such a large deviation cannot however be 
explained only by long range correlations which, as depicted in Fig. II. 6[ 
involve excitation energies of tenth of MeV. It has been argued that the 
deviation is evidence of SRC which, by populating highly excited states, 
strongly deplete shell model states below the Fermi sea. 
It should be pointed out that, for decades, understanding the role played 
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by SRC in nuclei has been a very elusive problem of nuclear physics, due 
to the difficulties in isolating the signal of SRC. As a matter of fact, at low 
and medium energies, effects due to the final state interaction of the knocked 
out nucleon with the residual system, as well as effects from the excitation 
energy of hit nucleon, could produce the same final state as the one which is 
expected to be produced by an initially correlated pair. Investigating SRC 
thus requires high energy probes, in order to cover kinematical regions with 

g2 > 1 [GeVjcf , xb, = > 1 (1-29) 

which are available in modern lepton and hadronic accelerators like, e.g., 
JLab (USA), GSl (Germany) and JPARC (Japan). In Eq. fOOj) . is the 
four- momentum squared of the virtual photon, u is the energy transfer, rriN 
is the nucleon mass, and Xsj is the Bjorken scaling variable. 
It is only recently that such experiments could be performed, and the exci- 
tation strength due to SRC could be observed (cf. Fig. II. 6p . 



1.2.2 Exclusive lepton scattering: the A{e,e'pN)X re- 
action 

Experimentally, high momentum probes can knock out a proton off a nucleus, 
leaving the rest of the system nearly unaffected. If, on the other hand, the 
proton being struck is part of a correlated pair, the high relative momentum 
in the pair leads the correlated nucleon to recoil and be ejected and detected, 
as pictorially shown in Fig. 11.81 This triple coincidence experiments have 
been performed on targets, both at Brookhaven National Laboratory 
(BNL) using incident protons [HI 1121113], and at Jefferson Laboratory (JLab) 
by using electrons [5]. 

The BNL experiment EVA (E850) studied the process 

p+^^C^2p + n + X (1.30) 

with proton beam momentum p ranging from 6 GeV/ c up to 9 GeV/ c. 

In the final state, the scattered proton, with momentum p', and the knocked 

out proton, with momentum 

P/^q = p-p' (1-31) 

where q is the momentum transfer, are detected. By assuming that the 
knocked out proton leaves the nucleus without interacting with the residual 
system {A — 1), one has 

p/ = p. + q- (1-32) 
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correlations in nuclei 




□itudrrof n^iiLron ^^^^^^^ 



Figure 1.8: A simple cartoon to illustrate the A{e,e'pN)X reaction. The 
incident electron couples to a nucleon-nucleon pair via a virtual photon. In 
the final state, the scattered electron and struck nucleon are detected along 
with the correlated nucleon that is ejected from the nucleus [S]. 



By introducing the so-called missing momentum 



q-p/ 



one has 



A-l 



-Pi 



as shown in Fig. 11.91 

If the knocked out nucleon was initially correlated with a neutron, 
nucleus {A — 2) almost at rest, one should observe and detect the 
neutron with momentum 



Pn 



-Pi = Pr 



Within these assumptions, by plotting the momentum p„ of the 
neutron as a function of the cosine of the angle between and p„. 



cos 7 



IPillPn 



(1.33) 

(1.34) 

with the 
recoiling 

(1.35) 

recoiling 
i.e. 

(1.36) 



15 



Chapter 1. Short range correlations in nuclei 




(b) 



Figure 1.9: Two simple cartoons to illustrate the ^^C{p,p'pn)X reaction stud- 
ied at BNL. (a) The incident and the struck protons, with initial momenta 
p and Pj, and final momenta p' and pj, respectively, are detected in coinci- 
dence with a neutron with momentum p„; (b) The same process, viewed also 
in terms of the missing momentum p^ and the three-momentum transfer q. 
The notation for angles is self explaining. 
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which is known as the number of directional correlations, one should observe, 
at high values of the momentum p„, a strong back-to-back directional corre- 
lation between Pj and p^, due to SRC, which has indeed been observed, as 
shown in Fig. 11.101 
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Figure 1.10: The correlation between p^ and its direction 7 relative to Pj. 
Data labeled by 94 and 98 are from Refs. [ITl WI\ . The dotted vertical line 
indicates the Fermi momentum kp = 0.221 GeV/c. After Ref. 



(1.37) 



In Ref. Sal m], the ratio 

Number of (p,2pn) events(pi,p„ > kp) 

r = . 

Number of (p,2p) events(pj > kp) 

has bee extracted from these data. Eq. fll.37p represents the measure of 
correlation of backward neutrons with initial momentum p„ ~ when 
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IPnMPil ^ ^f; by determining the ratio fl 1.371) it has been demonstrated 
that 921^8% of protons in with momenta pi > 275 MeV/c are partners 
in n-p SRC pairs. 



10 
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Figure 1.11: The fractions of SRC pair combination in ^^C, resulting from 
{e,e'pp) and {e,e'pn) reactions performed at JLab [5j, as well as from BNL 
{p, 2pn) data [HI S^l US]. The results show the dominance of p-n pairs [5]. 



The Jefferson Lab experiment has demonstrated that nearly all nucleons in 
^^C with momentum in the range 300-600 MeV/c have a correlated nucleon 
partner with roughly equal and opposite momentum ES]- By comparing 
n-p to p-p pairs yields, it has also been found that SRC are mainly due 
to n-p pairs, whose probability results of the order of 18 ± 5, as shown in 
Fig. 11.111 [5j. Calculations explain the magnitude of this n-p to p-p ratio 
as arising from the short range tensor part. Both experiments have shown 
that recoiling nucleons, with a momentum above the Fermi sea level in the 
nucleus, are part of a correlated pair, and both observed the same strength 
of p-n correlations. This confirms that the process is accessing a universal 
property of nuclei, unrelated to the probe [15]. 
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1.2.3 Inclusive lepton scattering: the ^4(6, e')X reaction 

The high energy A{e, e')X reaction depicted in Fig. I1.12[ i.e. the process in 
which only the scattered electron is detected, represents the simplest reaction 
to investigate SRC and, in particular, to measure the probabilities of SRC in 
nuclei. 




Figure 1.12: A simple picture of the A{e,e')X process. An electron with 
momentum k knocks out a nucleon N with final momentum pj and, after 
interaction, is detected with momentum k'. The knocked out nucleon is not 
detected. 

In Ref. [6|, [7] it has been supposed that, at high momentum transfer, the 
inclusive cross section off a nucleus A can be written as follows 

i=2 ^ 

where crA{xBj,Q'^) is the electron-nucleus cross section, written in terms of 
the 4- momentum transfer and the Bjorken scaling variable xbj, o'j{xBj,Q'^) 
is the inclusive cross section off a correlated cluster of j-particles, and, even- 
tually, the quantity aj{A) is the probability of finding a nucleon in the cluster 
j. By this way, one expects that the inclusive cross section is dominated, at 
1 ^ ^fij ~ 2, by the absorption of the virtual photon on a pair of corre- 
lated nucleons and by the elastic rescattering in the continuum, whereas, at 
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2 < < 3, it is governed by three-nucleon correlations (3NC), and so on. 
In the process under analysis, as shown in Fig. \1.12\ an electron beam with 
momentum k, knocks out from the nucleus A, a nucleon with final four- 
momentum given by the following energy and momentum conservation 

pj = {q + P^- PA-if = ml (1.39) 

where q is the four- momentum transfer, Pa and Pa-i are the momenta of 
the target nucleus and the residual system {A — 1), respectively, and m-Ar is 
the nucleon mass; as already pointed out, after interaction only the scattered 
electron, with momentum k', is detected. In Fig. I1.13[ the minimum value 
of the missing momentum 

Pm = q-P/ = PA-i (1.40) 

defined in terms of the the three-momentum transfer q = k — k' and the 
momentum of the residual system P^-i, is plotted versus the Bjorken scaling 
variable XBj- It can be clearly seen that, for any nucleus A and fixed Q^, 
the minimum recoil momentum |p™"| increases with XBj-, and exceeds the 
average Fermi momentum in nucleus A at xsj > the latter depending 
upon the nucleus [16]. As already pointed out, SRC correspond to high 
momentum components of the nuclear wave function, therefore, with the 
gradual increase of x^j, the virtual photon should first probe high momentum 
configurations due to 2NC and then, following Eq. fll.38p . at Xsj > 2, it 
should probe 3NC [6l[7]. 

To avoid the difficulties due to the theoretical calculations of electron off-shell 
nucleon cross section, inclusive data from JLab-HallB have been analyzed 
[H HB] not directly in terms of cross sections, but by taking the cross section 
ratios of '^He, ^^C and ^^Fe to ^He. 

The experimental cross section ratio of nucleus A to the nucleus ^iJe 

plotted versus the Bjorken scaling variable xbj, is shown in Fig. 11.141 In Eq. 
(I1.4ip . a{A) and ai^He) are the A{e, e')X and ^He{e, e')X inclusive cross sec- 
tions, respectively, and cXep and are the elementary elastic electron-proton 
and electron-neutron scattering cross sections. Three different kinematical 
regions are clearly seen: 

• XBj ^ 1-5: here the shape of the ratio is governed by the different 
behavior of the magnitude of the quasi elastic peak in different nuclei 
(higher peaks for light nuclei, and lower peaks for heavy nuclei); 
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1 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

Figure 1.13: The minimum recoil momentum versus the Bjorken scaling vari- 
able XBj for a) Deuteron, calculated at several [Gel^^],and for b) different 
nuclei at = 2.0 GeV^. Horizontal lines at 250 MeV/c indicate the Fermi 
momentum typical of uncorrelated motion of nucleons in nuclei ^6] . 

• 1-5 < a;Bj < 2: this plateaux is interpreted as evidence of two-nucleon 
correlations (2NC), which in complex nuclei and in ^He should differ 
only by a scale factor; 

• 2 < Xfij < 3: the presence of a second plateaux is interpreted as 
evidence of 3NC. 

The relative per-nucleon probabilities a2{A,^ He) and a^lA,^ He) of finding, 
respectively, 2NC and 3NC in nuclei relative to ^He, extracted from the 
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Figure 1.14: Weighted cross section ratios of a) ^He, b) ^^C and c) ^^Fe to 
^He as a function of xbj for > 1.4 GeV^. The horizontal dashed lines 
indicate the 2N and 3N scaling regions used to calculate the per-nucleon 
probabilities of 2N and 3N SRC in nucleus A relative to ^He. After Ref. 
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experimental results shown in Fig. 11.14^ are listed in Table 11.11 Here also 
shown are the absolute (per-nucleon) values a2{A) and as (A) of the same 
probabilities in nucleus A, calculated in Ref. jl] from 



= (1-42) 



«3(A^^e) = ^^^ (1.43) 
by using realistic wave functions of ^He and Deuteron nuclei. 



A 


a2{A,^He) 


a2N{A) (%) 


a3{A,^He) 


a3N{A) (%) 


3 


1 


8.0 ± 1.6 


1 


0.18 ±0.06 


4 


1.93 ±0.01 ±0.03 


15.4 ±3.2 


2.33 ±0.12 ±0.04 


0.42 ±0.14 


12 


2.49 ±0.01 ±0.15 


19.8 ±4.4 


3.18 ±0.14 ±0.19 


0.56 ±0.21 


56 


2.98 ±0.01 ±0.18 


23.9 ±5.3 


4.63 ±0.19 ±0.27 


0.83 ±0.27 



Table 1.1: The relative per-nucleon probabilities a2{A,'^ He) and a^lA,"^ He) 
of 2NC and 3NC in nucleus A relative to ^He, and the absolute value a2N{A) 
and a3N{A) of the same probabilities in nucleus A (in %), from Ref. 
Errors shown are statistical and systematic for 02 and 03, and are combined 
(but systematic dominated) for a2N and a^j^. 



It should be pointed out that no direct calculations of the inclusive cross 
section ratio shown in Fig. 11.141 have been performed so far. These calcu- 
lations would represent a relevant contribution towards the solution of the 
longstanding problem concerning the role played by SRC in nuclei. 
In this Thesis, preliminary results of the calculation of the inclusive ratio 
ri^^Fe/^He) will be given in Chapter 13 and a new approach [17] to inclusive 
electron scattering at high momentum transfer will be illustrated. 



1.3 Relevance of short range correlations in 
various fields 

As shown in Fig. II. 5[ SRC can lead, in small portion of time, to local den- 
sities in nuclei comparable with those in the core of neutron stars, i.e. up to 
~ 7 times the average nuclear density po ~ 0.17 N/ fm^ . Moreover, combin- 
ing the results of the various experiments we have illustrated, the probability 
of finding nucleons in ^^C which move in an average potential has been found 
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IZ] Single nucfeons 

n-p I n-n CH p-p 



Figure 1.15: Probability of finding shell model and correlated nucleons in 
After Ref. [5]. 



to be of the order of ~ 80%, whose ~ 60-70% due to shell model potential, 
whereas ~ 10-20% due to long range correlations; the remaining 20% repre- 
sents the SRC contributions, which are dominated by p-n correlations, due 
to the strong tensor force which acts between a n-p pair and does not act in 
n-n and p-p pairs, as shown in Fig. 11.151 [5] . 

Obtaining information on SRC phenomena in nuclei would have a strong im- 
pact on various fields of physics, e.g. in particle-, nuclear- and astro-physics. 
Let us briefly illustrate some of them: 



- NN interaction : NN interaction processes at large and intermediate 
distances are well described in terms of meson exchange, as shown in 
Fig. I1.16[ but such a picture does not allow a proper description of 
the strongly repulsive NN interaction at short distances. Therefore 
understanding SRC in nuclei should allow a deeper knowledge of the 
NN force gH]. 
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- Cold dense nuclear matter : the density between correlated nucleons is 
similar to the one in cold dense nuclear objects, as neutron stars, thus 
a deep understanding of SRC effects should lead to the formulation of 
a realistic equation of state for such systems. As shown in Fig. I1.17[ 
the core of neutron stars, neglecting SRC, could be well approximated 
by two independent Fermi gases: the prevalent one constituted by neu- 
trons, and the smallest one by protons. Even if the number of protons 
is small, the strong n-p correlation acts towards the coupling of the 
two Fermi liquids, thus affecting the equation of state. Moreover, it 
becomes worth analyzing the role played by SRC phenomena in the 
URCA processes 

n ^ p + e + Ue (1-44) 
e + p-^n + Ue (1-45) 

which involve neutrino and antineutrino emission, leading to 
changes in the physical properties of the system [2]. 
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Figure 1.17: Illustration in momentum space of the momentum distribution, 
kpi oi protons and neutrons in neutron stars. The left side shows the case 
where the protons and neutrons weakly interact and can be approximated as 
separate Fermi spheres, with the neutron momentum much greater than the 
proton momentum. The figure on the right shows how strong neutron proton 
SRC cause protons to escape their Fermi sphere and have higher momentum 
then would otherwise be allowed. After Ref. [i5] . 

- Hadron properties: another important problem in modern nuclear physics 
is given by the modifications of hadron properties such as masses and 
radii in the nuclear medium. The nuclear EMC effect, i.e. the change 
in the inclusive deep inelastic structure function of a nucleus relative 
to that of a free nucleon, induces such modifications, which are of fun- 
damental importance in understanding implications of quantum chro- 
modynamics (QCD) for nuclear physics. As recently shown, possible 
modifications of nucleon properties induced by the medium can be bet- 
ter studied by analyzing the short range properties of nuclei |19] . 

- QCD: another open problem in modern nuclear and particle physics 
which could be answered by investigating SRC, concerns the role played 
by quark and gluon degrees of freedom in nuclear matter at distances, 
as depicted in Fig. I1.18[ which are expected to be most relevant in this 
interaction range. 

- High energy scattering processes: recently it has been shown [501 EI] 
that high energy diffraction effects in nuclei are strongly affected by 
SRC, which appear to be of the same order as Gribov inelastic shadow- 
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Figure 1.18: Distances between correlated nucleons are so close that non 
nucleonic degrees of freedom should reveal themselves. 

ing [52]. An example is given in Fig. ll.lQI by the total neutron-Nucleus 
cross section. 

A deep comprehension of SRC phenomena is mainly addressed to answer the 
following questions: 

• What is the percentage of finding correlated nucleons in nuclei? 

• What is the relevance of three nucleon SRC in nuclei? 

• What is the ratio of pp to nn pairs? 

• Are tensor forces relevant for SRC? 

• Are these nucleons different from free nucleons? 

• Which type of analysis should be used in order to obtain information 
on SRC? 

The increasing interest of the scientific community towards the comprehen- 
sion of SRC in nuclei can be stressed by the large number of conferences 
and workshops organized in the last year. Some examples are shown in Figs. 
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Figure 1.19: The neutron-nucleus total cross section, for ^ife, ^^C, ^^O and 
^°^P6. Lejt 'panel: the result without the inclusion of correlations; dotted 
curves: shell model contribution; dot-dashed curves: shell model contribution 
plus inelastic shadowing effects. Right panel: results with the inclusion of 
correlation. Dotted curves: shell model contribution; dashed curves: shell 
model contribution plus two-nucleon correlations; solid curves: shell model 
contribution plus two-nucleon correlations plus inelastic shadowing effects. 
Experimental data from [53l|5l]. After Ref. |50] . 
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Figure 1.20: 2009 Mini-symposium on Nuclear Structure at Short 
Distances I, held within the American Physical Society April Meeting, May 
2 - 5, 2009, Denver (USA) [IZl [55]. 
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KEK theory center workshop on 

Short-range correlations and tensor structure at J-PAEC 
Sept. 26, 2003 (Friday), KEK, Building 4, Eoom 3AM-6PM 
UEL: http://www-oonf.l[ek.jp/hadronl/j -paro-sro09/ 

Orgainsers: Shunso Kumano (KEK), Hike Merita (Sapporo Gakuin Univ), 
Shin'ya Sawada (KEK) 



Program 



8:30 - 3:00 Eegistration 

Chair: Osamu Morimatsu (KEK) 

3:00 - 3:55 (50+5) Claudio Ciofi degli Atti (Perugia Univ) 

Short Eange Correlations and their impact on nuclear and particle physics 
and astrophysics: recent advances and possible studies at J-PAEC (pdf ) 

3:55 - 10:25 (26+5) Shin'ya Sawada (KEK) 

Status of J-PAEC facility and hadron hall (pdf) 



10:25 - 10:50 treak 



Chair: Koiohi Saito (Tokyo Univ. Science) 

10:50 - 11:10 (16+5) Chiara Benedetta Messetti (Perugia Univ) 

Two- and three-nuoleon correlations and inclusive electron scattering (ppt) 

ll:lO - 11:35 (20+5) Hiko Merita (Sapporo Gakuin Univ) 

Nulcear spectral function based on the two-nucleon correlation model (ppt) 

Figure 1.21: Short-range correlation and tensor structure at J- 
PARC, September 25, 2009, KEK Tsukuba (Japan) [56j. 
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Figure 1.22: The Jefferson Laboratory Upgrade to 12 GeV, September 
14 - October 16 & October 26 - November 20, 2009, INT Seattle (USA) [5?]. 
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Figure 1.23: Probing Hadron Structure from Hard Exclusive Pro- 
cesses, November 9 — 11, 2009, Garching (Germany) |58j . 



32 



Chapter 2 



The spectral function and the 
nucleon momentum 
distributions: two- and 
three-nucleon correlations 



Introduction 

In order to introduce our new approach to the inclusive cross section, it 
is necessary to recall some basic concepts and general properties about the 
spectral function and the nucleon momentum distributions. 

2.1 The spectral function 

The one-body spectral function, which is defined as follows (see e.g. [59] ) 
P^'ik^E) = < "^^lal 5 {E-{H-EA))a^^\^ A > (2.1) 



represents the joint probability distribution of finding in the target nucleus 
a nucleon with momentum A; = |k| and removal energy 

E = E^i^ + E\_^ = mN + Ma-i - Ma + E\_^ . (2.2) 

Here, aj^^ and creation and annihilation operators of a nucleon with 

momentum k and spin a, respectively; H is the intrinsic Hamiltonian of A 
interacting nucleons; is the ground state eigenfunction with eigenvalue 
Ea, total angular momentum J and projection M; eventually, rriN is the 
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mass of the nucleon, M4 the mass of the target nucleus, and M^-i the mass 
of the system (A — l), with excitation energy In what follows, for ease 

of presentation, the absolute value of a vector a will be indicated as a = |a|. 
By placing in Eq. (12.11) the completeness relation 



Ei^A-i><^Lii = i 

the spectral function becomes 



(2.3) 



M,a f 



(27r)3 (2J 
X 6(^E-{E{_,-Ea)) 



(2.4) 



where is the intrinsic eigenfunction of the final state / of the Hamilto- 

nian Ha-i with eigenvalue 



A-l\ 



E 



A-l 



and Gf {z) is the overlap integral 



Gf^iz) =< xV\ . . . , y)|^A(x, . . . , y, z) > 



(2.5) 



(2.6) 



1/2 

where is the two-component Pauli spinor of the nucleon. Since the set 
of states / also includes the continuum states of the residual (A — l)-nucleon 
system, the sum over / in Eq. (12. 4p stands for summation over the discrete 
states of the (A — l) system and integration over the continuum states. Thus 
the spectral function exactly includes all final states interactions in the states 
of the (A — 1) system, the only plane wave being that describing the relative 
motion of the knocked out nucleon and the (A — l) system. 
The spectral function obeys the normalization condition 



P^{k,E)dk dE 



(2.7) 



and, owing to Eq. (12. ip . it can be written as follows 

P^(A;,E) = P(f(A;,E) + Pi^(A;,E) 
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where the contributions from different final nuclear states have been explicitly 
separated out, namely 



^ I ^ 'I M,<T f<c ^ 

X S(^E- {e{_^ - Ea)) (2.9) 
includes the ground and one-hole states of the {A — l)-nucleon system, and 

(2.)3(2J + 1) J 



X S [e - {Ei_, - Ea)) 



(2.10) 



more complex highly excited configurations generated in the target ground 
state by NN correlations. Here / < c (/ > c) means that all final states of 
the residual system below (above) the continuum threshold are considered. 



2.2 The nucleon momentum distribution 

By definition, the calculation of the spectral function requires the knowl- 
edge of the whole set of final states "^a-i- The calculation of the nucleon 
momentum distribution requires, on the contrary, only the knowledge of the 
ground state wave function. As a matter of fact, the momentum distribution 
is defined as follows 

n^{k) j dz dz' e*'^-^"-"') p(z, z') (2.11) 

where 

p(z,z') = Jd^...dy [vl/^(x...y,z)]*vl/0(x...y,z') (2.12) 
is the non diagonal one-body density matrix. 

The spectral function and the momentum distributions are related by the 
momentum sum rule 

r+oo 

n^{k)^47r / P^{k,E)dE (2.13) 

J E 
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which is readily obtained by inserting the completeness relation f l2.3p in Eq. 
Within the decomposition rule given by Eq. (12. Sp . one has 



(2.14) 



with 



ntik) = 4vr / P,^{k, E)dE = — — — - / e'^^ (z)ciz 

J Emin V 'Ma 



(2.15) 



and 



f+OO -1-1 /■ 

ntik) = 4vr / P,^{k, E)dE=— ——- E / ^^"'^ Gf-iz)dz 



(2.16) 

The integral of the nucleon momentum distributions yields the spectroscopic 
factors (or occupation probabilities) 



and 



5*0 



Si = 



dk k"^ n^ik) 



dk n^{k) 



which, owing to the normalization condition 

/ dk P n^{k) = 1 
Jo 

satisfy the relation 

So + Si = l. 



(2.17) 



(2.18) 



(2.19) 



(2.20) 



Note that, in the independent particle shell model description, as will be 
discussed in more detail in §2.31 the spectral function P^{k, E) can be written 
as follows 

(2.21) 



where 



/• + 00 

/ k"^ dk na{k) = 1 
Jo 



(2.22) 



and Aa is the number of nucleons in the state a {A = with removal 

energy and nucleon momentum distribution Uaik). 
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10 
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1 



3 



k [fm '] 



Figure 2.1: The nucleon momentum distribution n{k) = n^^{k) vs. the 
momentum k. The dotted line is the mean field contribution, the dashed line 
the one arising from central forces, and the long line includes both central 
and tensor forces. The full dotts are the results of the Variational Monte 
Carlo approach [33]. After Ref. [31] . 

In Fig. 12. 1^ the correlated nucleon momentum distributions resulting from 
the cluster expansion techniques discussed in §1.11 [31] are compared with the 
mean field component and the variational Monte Carlo calculations [33], in 
case of ^^O. It can be clearly seen that the mean field distributions almost 
totally exhaust the low momentum part of n"^{k), and drop to zero at k > 
1.5— 2/m~^; on the contrary, because of NN correlations, the high momentum 
tails are entirely governed by tensor forces, acting in T = and S = 1 states 
which, at high momenta, are several orders of magnitude larger than the 
predictions from shell model calculations. Moreover, as shown in Fig. 
12.21 and, more clearly, in Fig. 12. 3[ apart from a scaling factor, the nucleon 
momentum distributions n-^^k) seem to be almost independent of the atomic 
weight A] thus, a simple model for the nucleon momentum distributions 
n^{k) could be written as follows [60j 




(2.23) 



(2.24) 
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Figure 2.2: The nucleon momentum distributions n^(A;) for nuclei ranging 
from '^H to Nuclear Matter {NM). It can be seen that, at high values of the 
momentum A;, n'^ik) can be considered as a rescaled version of the momentum 
distributions of "^H. After Ref. [BUI E5]. 



where, because of the effects of correlations, which depopulate states below 
the Fermi sea, one has 

/•+00 

/ Pdkn^{k)<l (2.25) 
Jo 

and 

r+oo 

/ k^dknf{k)>0. (2.26) 
Jo 

Here the low momentum component, owing to f l2.13p . is nothing but Eq. 
( I2.2ip integrated over the removal energy E, with the shell model momen- 
tum distribution n-^ replaced by the modified ones which takes into account 
the effects of correlations; on the contrary, the high momentum part is noth- 
ing but the Deuteron nucleon momentum distribution n^{k) rescaled by a 
constant C^, which depends on the nucleus A under consideration, and whose 
values are listed in Table 12.11 We call such a behavior Deuteron scaling. 
Deuteron scaling is even more pronounced in the the two-nucleon momentum 
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Figure 2.3: The ratio n^{k) to the Deuteron momentum distribution n^{k). 
After Ref. (601. 



Nucleus 




4//e 






40Ca 


56pg 


208p^ 


iVM 




1.9 


3.8 


4.0 


4.2 


4.4 


4.5 


4.8 


4.9 



Table 2.1: Values of the constant appearing in Eq. (1^ 



distributions, defined as follows 



n^i^Hki,k2) = J rfziciz2ciz;dz^e' ^^(^^-<)e*''^(^^-^'2)p^^^^ (zi, Z2, z;, z'^) 

(2.27) 

which have recently been calculated in Refs. [9l l3T] . 

The knowledge of (ki, k2) allows one to calculate the relative and center 
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of mass (CM) momentum distributions of a two- nucleon pair A^i N2-, i.e. 

n^ei (k.ez) = j d\.cM (k.,, + -k.,, + ^) (2.28) 

n^c'M^ (kcM) = j dKei (k.,, + -k.,, + ^) (2.29) 

where 

krel = '~ ' kcAf = ki + k2 (2.30) 

are the relative and CM momenta of the pair N1-N2, respectively, with ki 
and k2 being measured from the CM of the system, and 

p^^^Mzi,Z2,z;,z'2)^ I rfx...y [vI/° (x...y,Zi,Z2)]* n(x...y,z;,z'2) 

(2.31) 

is the off-diagonal two-body density matrix. 

Deuteron scaling can be clearly seen in Figs. 12. 4^ 12.51 and 12. 6^ where the 
quantity n(krei, \^cm = 0), describing back-to-back nucleons, is plotted versus 
the relative momentum hrei- These figures clearly illustrate the dominant 
role of the tensor forces in producing a substantial difference between p-n 
and p-p two-nucleon momentum distributions, both in few-nucleon systems 
and complex nuclei. 
In Ref. [To], the following ratio 



/a dkrelk'^gi [npp{krel-, 0) + Upnikreh 0)] 

which represents the percentage probability of finding a p-n pair in the nu- 
cleus A, has been calculated. The results, listed in Table [2^ show that PpN 
is proportional to the percentage of p-N pairs, when integration runs over 
the whole range of kreu on the contrary, when integration is limited to the 
correlation region, the percentage of p-n pairs results much larger than that 
of p-p pairs, which is a clear consequence of the effects of the tensor forces 
acting between a proton and a neutron. 



2.2.1 The saturation of the momentum sum rule 

A relevant relationship between high momentum and high removal energy 
components can be obtained by considering the partial momentum distribu- 
tion t62j 

nj{k)=47i / dEP^{k,E) (2.33) 

J E 
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1 2 3 4 5 

Figure 2.4: The n-p (hnes) and p-p (symbols) two-nucleon momentum distri- 
butions Pnn{(1, Q) = n{kreh kcKi) in various nuclei as functions of the relative 
momentum q = k^ei at vanishing total pair momentum Q = kcM- After Ref. 
0. 
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Figure 2.5: The n-p (lines) and p-p (symbols) momentum distributions in ^He 
obtained with different Hamiltonians. Also shown is the scaled momentum 
distribution for the AV18 Deuteron; its separate S- and D-wave components 
are shown by dotted lines. After Ref. [9]. 
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Figure 2.6: The two-nucleon momentum distributions n{krei, kcu = 0) vs. 
the momentum k^eu for (upper panel) and ^°Ca (lower panel). After 
Ref. [10]. 
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A 


Ppp (%) [0, oo] 


Pp„ (%) [0, cx)] 


Ppp (%) [1.5, 3.0] 


Ppn (%) [1.5,3] 


4 


19.7 


81.3 


2.9 


97.1 


12 


30.6 


69.4 


13.3 


86.7 


16 


29.5 


70.5 


10.8 


89.2 


40 


31.0 


69.0 


24.0 


76.0 



Table 2.2: The p-p and p-n percentage probability given by Eq. fl2.32p eval- 
uated in the momentum range shown in square brackets in fm~^. After Ref. 



where the upper limit of integration Ef can be varied from Emm to oo. By 
definition, Eq. fl2.33p represents that part of n^{k) which is due to final 
{A — l)-nucleon states with E < Ef. In the limit Ef ^ oo one gets 

nj{k) n^{k) (2.34) 

and the momentum sum rule given by Eq. (I2.13P is recovered. Thus the 
behavior of n^{k) as a function of Ef provides information on the saturation 
of the momentum sum rule and the relevance of binding effects. 
The saturation of the momentum sum rule for ^He [62] and Nuclear Matter 
[65] is shown in Fig. 12. 7[ using realistic spectral functions. It can be clearly 
seen that, at < 1.5 the momentum sum rule is saturated already at 

values of Ef very close to -Emin, whereas at A; > 1.5/m~\ on the contrary, the 
momentum sum rule is saturated only when high values of Ef are considered. 
This can be explained in terms of the spectral function P^{k, E) appearing in 
Eq. fl2.33p which, at low momenta, is dominated by its component P^ik, E), 
whose strength is almost totally concentrated at low values of the removal 
energy; at high momenta, on the contrary, it depends upon P^{k, E), spread, 
for a given k, over a wide range of values of E. 



2.2.2 Probabilities of independent particle and corre- 
lated momentum components 



Let us define 

"+00 



So= dk^n^ik) (2.35) 

Jk* 



and 

"+00 



Si = I dkk^n^ik) (2.36) 
Ik" 
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12 3 4 

k (fm'^ ) 




12 3 4 

k (fm''' } 

Figure 2.7: The saturation of the momentum sum rule in '^He (a) and infi- 
nite nuclear matter (b). The dotted and solid lines correspond to the mo- 
mentum distribution nQ(k) and to the total momentum distribution n^(k), 
respectively. In case of ^He the dot-dashed, dashed, and long dashed lines 
correspond to Eq. (12.331) calculated at Ef = 17.75, 55.5, 305.5 MeV^, whereas 
for nuclear matter the dot-dashed and dashed lines correspond to Ef = 100 
and 300 MeV, respectively. The spectral function for ^He is from Ref. [62j . 
and for nuclear matter from Ref. [63[ IM] . After Ref. [60] . 

as the probability of finding independent particle and correlated momentum 
components in nucleon momentum distributions, with k* ranging from to 
oo. When A;* = 0, the spectroscopic factors given by Eqs. (12.171) and (I2.18P 
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are recovered. 

The values of these two quantities, calculated for different values of k* and 
for different nuclei, are listed in Table [2731 It can be seen that: 







12^ 




k* [fm-1] 


So 


Si 


So 


^1 


So 


Si 


0.00 


0.80 


0.20 


0.80 


0.20 


0.80 


0.20 


0.25 


0.75 


0.19 


0.78 


0.19 


0.78 


0.19 


0.50 


0.55 


0.17 


0.69 


0.18 


0.69 


0.19 


0.75 


0.31 


0.14 


0.48 


0.15 


0.45 


0.17 


1.00 


0.14 


0.12 


0.24 


0.13 


0.17 


0.15 


1.50 


1.5 10-2 


8.3 10-2 


2.1 10-2 


8.510-2 


2.710-3 


0.11 


2.00 


7.4 10-^ 


6.1 10-2 


6.2 10-^ 


6.1 10-2 


2.5 10-6 


7.5 10-2 



Table 2.3: Eqs. fl2.35p and fl2.36p calculated at different values of the mo- 
mentum k*, for ^He, ^"^C and ^^Fe. 



• when the whole range of the momentum k is considered, the probability 
of finding independent particle low momentum shell model components 
is much larger than the probability to find correlated high momentum 
components; 

• for 0.75 ^ k* ^ 1.50 fm~^ the contribution from shell model and 
correlated nucleons is comparable; 

• for k* > 1.50/m-^ the probability of finding shell model nucleons is still 
different from zero, but orders of magnitude less than the probability 
due to correlated nucleons. 

Let us now discuss in more detail the low and high momentum components 
of the nuclear wave function. 



2.3 Low momentum components and the mean 
field structure of the nuclear wave func- 
tion 

Eq. (12. 9 p yields the probability distribution that the final {A — l)-nucleon 
system is left into its ground state, i.e. with Ej^_^ = and E = Emin- 
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As already pointed out, the shell model spectral function can be written as 
follows 

PsMik,E) = ^5^A„nf ^(A:)5(E- |ej) (2.37) 

a 

with Aa denoting the number of nucleons in the state a with removal energy 
Sa and nucleon momentum distribution na{k). At the same time, one has 

Pf{k,E) = (2.38) 

In Eq. fl2.37p . the sum over a runs only over the hole states of the target, 
which means that the occupation probability Sa is 

The main effect of NN SRC is to deplete states below the Fermi level and 
to make the states above the Fermi level partially occupied. By such a 
mechanism 

Po^(fc,i?)^P#M(fc,^) (2.40) 

and 

P^ik,E)^0 (2.41) 

Disregarding the finite width of the states below the Fermi level, the modified 
shell model contribution can be written as 

PsMik, E) = ^J2 ^"^"(^) - (2.42) 
where the occupation probability for hole states is 

POO 

S^= I k'^dkn^ik) < 1 (2.43) 
Jo 

Eq. 02.421) drops down very quickly for k > kp and E > Ep, where the 
spectral function behavior is governed by NN correlations, and thus by its 
correlated component P^{k,E). 



2.4 Two-nucleon correlations 

In Fig 



2.8 a 



the spectral function of ^He corresponding to the AV18 inter- 
action [20 j is shown versus the momentum k and the removal energy E. A 
typical feature of the spectral function, common also to the spectral function 
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Figure 2.8: (a) The ^He spectral function corresponding to the AV18 poten- 
tial [20]; (b) Some cuts of the spectral function P{k,E) at high values of k 
and E. After Ref. [MIEZ]. 



of Nuclear Matter [63| IM] . is that it exhibits, at high values of k and E, 
broad peaks located at 

E*A 1 - ^7-^ ^ (2.44) 

A- I 2mjv 



and whose width increases with /c, as clearly exhibited in Fig. 2.8(b) , where 
some cuts of P^{k, E) are shown. Eq. fl2.44p is generated by 2NC in nu- 
clei, as will be demonstrated in what follows, by illustrating different model 
developed in the last few years, differing in the description of the nucleon 
configurations inside the nucleus. To this end, in Fig. 12. 9[ a simple cartoon 
of a target nucleus A and its nucleon constituents is given. 

2.4.1 The naive two- nucleon correlation model 

In Ref. [HI [7], a first microscopic model leading to Eq. fl2.44p has been 
proposed. It is based upon the assumption that, referring to Fig. I2.10[ the 
high momentum ki = k of nucleon '1' is entirely balanced by the momentum 
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Figure 2.9: A momentum description of a nucleus A. Nucleon '1', also called 
the "active" nucleon, has momentum ki = k, and the nucleus {A — 1) consists 
of nucleon '2', with momentum k2 and nucleus {A — 2), with momentum 
Ka_2; the latter is defined in terms of nucleon '3' and nucleus {A — 3), with 
momenta ks and K^_3, respectively. 



k2 ^ — k of nucleon '2', whereas the residual spectator system has total 
momentum Ka-2 — 0. Energy conservation for such a mechanism requires 
that 



where 



1 ~ — r-^^-T (2.46) 



is the recoil energy of the residual {A — l)-nucleon system, whose intrinsic 
excitation energy could therefore be written as 

E*A 1 - (2.47) 
A-l 2mN 
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Figure 2.10: The naive two-nucleon correlation model: the high momentum 
ki = k of nucleon '1' is entirely balanced by the momentum k2 ~ — k of 
nucleon '2', with the system {A — 2) being at rest, i.e. with momentum 
Ka-2 ^ 0. 



Within such a picture, the nucleon spectral function simply reads as follows 
P^M^'^^E-ESl-^^^) (2.48) 

where 

E[ll ^ \Ea\ - \Ea-2\ (2.49) 

is the two-nucleon breakup threshold, being \Ea\ and |-Eyi_2| the (positive) 
ground state energies of nucleus A and {A — 2). 



2.4.2 The convolution formula 

The naive 2NC model has been implemented in [601 EH], by assuming that 
the momentum ki = k of nucleon '1' is not fully balanced by the momentum 
k2 of nucleon '2', but also by the residual {A — 2) system, as depicted in 
Fig. I2.11j in other words, it is assumed that the spectator system {A — 2) 
moves with small momentum K.a-2 7^ 0. We will call such a configuration 
few nucleon correlation (FNC) configuration. 

Introducing the previously defined center of mass, kcM, and relative, krei, 
momenta of the correlated pair as follows 



kcA/ 

krei 



ki + ka : 
ki - ka 



k + ks 
k-k2 



(2.50) 
(2.51) 
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Figure 2.11: Tlie few nucleon correlation model: the high momentum ki = k 
of nucleon '1' is entirely balanced by the momentum k2 of nucleon '2', and 
by the momentum K.a-2 of the residual system {A — 2), with Ka-2 ^ |k2|. 



momentum conservation yields 



ki + ka + Ka-2 = 

Ka-2 = — kcM 



and energy conservation reads as follows 

E-E{_, + Ea = E + \Ea-2\ 
where 



If |2 

\^2,{A-2)\ 

2/i2,(A-2) 



\E, 



E-lE 



.(2) 
^thr 



(2.52) 
(2.53) 



If |2 

\^2,{A-2)\ 

2/i2,(A-2) 

(2.54) 



If |2 

1^2,(^-2)1 



(A -2) 



2Ai2,(A-2) 2mAr(A - 1) 



{A - 2)k2 - Ka-2 
(A-l) 



(2.55) 



is the energy of the relative motion of particles '2' and \A — 2)\ i.e. the exci- 
tation energy E\_^ of the residual {A — 1) system. In Eq. (I2.55p . eventually, 



f^2,{A-2) 



A-2 

A-l 



rriN 



(2.56) 



is the reduced mass of particles '2' and '(A — 2)'. 

Given these assumptions, the spectral function, at high k and high E, has 
been obtained in Ref. |60] by assuming that the ground state wave function 
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of nucleus A factorizes as follows 



(2.57) 

where 

X = ri - r2 y = rs - ^^—^ (^-^S) 

are the relative and CM coordinates, respectively, with = Ra-2, A is 
a proper antisymmetrization operator, and is a short-hand notation for 
the standard Clebsh-Gordan coupling of orbital and spin angular momenta; 
{$„(x)} and {Xm(y)} represent a complete set of states describing the rel- 
ative and CM motion of the pair, and, eventually, \E'^'l"2^({rj}A-2) is the 
complete set of states describing the {A — 2)-nucleon system. 
As already explained, SRC correspond to the high momentum and high re- 
moval energy components of the spectral function P'^iJv^E). In momentum 
space, SRC are described by a correlated pair with a very high relative mo- 
mentum krei > kp, and a low CM momentum kcM ^ kp; the latter condition 
justifies the choice of the CM motion in s-state, so that 



Ki{^^}A) ^ A I Xo(y) Yl ® *A-2({r.}A-2)] \ (2.59) 

njA-2 J 

where Xo(y) describes the low momentum CM wave function in s-state. Since 
the s-wave motion also implies that the system {A — 2) is in the ground state 
or in a low energy excited state, one can write, eventually. 



where 



Ki{^^}A) ^ A |xo(y) [$(X) ® K-2i{r^}A-2)] } (2.60) 

$(x) = J2 anoo'fn(x) (2.61) 



describes the relative motion of the correlated pair in the nuclear medium. In 
Eq. (I2.60p . fA-2 = denotes the excitation spectrum of the system {A — 2) 
which, since, as already pointed out, the CM of the pair involves only low 
momentum components, have been mainly limited to the ground state and 
to the (low-lying) excited states corresponding to configurations generated 
by the removal of two particles from different shell model states of the target. 
The above formalism leads to the following convolution formula, involving 
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only the relative and CM momentum distributions of the correlated pair 

^ ~ ^'^^ " 2M(A - 1) V (A^) 
X 5(k + k2 + KA-2) (2.62) 

where rirei and r?-^// are the momentum distributions of the relative and 
center of mass motion of the two nucleons in a correlated pair, respectively, 
and the other notations, following Fig. 12. 9^ are self explaining; eventually, 
the factor A/'2 satisfies the normalization condition given by Eq. (12. 7p . We 
stress once again that in such an approach, the residual {A — 2)-nucleon 
system is assumed to be in its ground state, which implies that the CM of the 
correlated pair moves with low momentum, so that only " soft" components of 
ncM contribute to the spectral function Ppj^c^ since we have seen in Chapter 
2 that high momentum components of a nucleon are linked to high excitation 
energies of {A — 1) and this has to be true also for the system {l-2)-{A — 2). 
The naive 2NC model is recovered by placing 

ncAii^CM) = 5(kcA/) (2.63) 

i.e. by assuming that the spectator nucleus is at rest. 

Integrating Eq. fl2.62p over k2 and the angular variables of K^_2, one gets 



P^j^cik, E) = M2 / dKA-2 Ka-2 UreliK) ncM{KA-2) (2.64) 



where 



= ^ \k±ko\ (2.65) 



ko = J2mr,^-^[E-Elll] (2.66) 



Ak^ + (A-2)B AKl 2 
= V 2(A-1) -4(J^ P-^^) 

The convolution formula fl2.64p has been calculated in Ref. ^0] by using the 
following effective relative and CM nucleon momentum distributions 

nJJ{krei) = C^nD{Kei) (2.68) 
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Figure 2.12: The nucleon spectral function of ^He \l2l and nuclear matter 
[551 El] versus the removal energy E for various values of the momentum 
k. For ^He (a) the squares, full dots, and open dots correspond to k = 
2.2, 2.8, 3.5 /m~^, respectively. For nuclear matter (b) the open squares, full 
dots, open dots, and full squares correspond to A; = 1.5, 2.2, 3.0, 3.5 /m~^, 
respectively. The solid line correspond to the FNC theoretical calculation, 
Eq. f ETMj) . After Ref. [SD]. 
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Figure 2.13: Tlie distribution of the cosine of the opening angle between 
the missing momentum and p^g^ = pcM — Pm? for the Pm = 0.55 GeV/c 
kinematics. The histogram shows the distribution of random events. The 
curve is a simulation of the scattering off a moving pair with a width of 
0.136 GeV/c for the pair CM momentum. After Ref. \38\ . 



respectively, with the parameter acM determined as explained in Ref. pU] . 
Recently, Urei and ucm have been obtained from many-body approaches 
in [To], separately for p-n and p-n pairs, and they quantitatively agree with 
the forms given by Eqs. fl2.68p and f l2.69p . In such a way, the convolution 
formula fl2.64p is completely defined in terms of many-body quantities. 
Eventually, it should be pointed out that, in the recent BNL experiment on 
^^C target [38], the CM momentum distribution has been measured, finding 
that it fully agrees with the prediction of Ref. [60] , as shown by the full line 



and 




(2.69) 



in Fig EH 
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2.5 Brueckner-Bethe-Goldstone theory and the 
validation of the convolution formula 

The convolution model can be microscopically derived from the Bruckner- 
Bethe-Goldstone (BEG) theory ^69j. Let us discuss this aspect in more detail. 
In nuclear matter the spectral function corresponding to the nucleon self- 
energy 

M{k,E) = V{k,E) +iW{k,E) (2.70) 
is given by the well known result 

jNM 



1 c^nn m 1 W{k,E) 



p''''\k,E) = — '^g{k,E) = - 2 

TT n [-E -ky2mN -V{k,E)f + W^{k,E) 

(2.71) 

where Q{k,E) is the single-particle Green function 

" -E - ky2mj, - V{k, E) - iW{k, E) ' ^^''^^^ 

It has be noticed that the real, V{k,E), and imaginary parts, W{k,E), of 
the self-energy are highly off-shell in the considered energy and momentum 
ranges. We are interested in the region where E is much greater than the 
Fermi energy Ep\ for high k and E, one has 

E + ^:^\V{k,E)UW{k,E)\ (2.73) 

and Eq. f l2.7ip becomes 

(2.74) 

where p denotes a "particle state" (outside the Fermi sea), and h{h') a "hole 
state" (inside the Fermi sea), with energies e(p), e{h) and e{h') respectively, 
and G is the BEG scattering matrix 

6*12(0;) = v + v ^ G{uj) = v + v- ^ 



U — Hq + tf] UJ — Hq + IT] 

+ V $ V § v + ... (2.75) 

00 — Ho + IT] LO — Hq + 17] 

where Q is the Pauli projection operator, which restricts the two nucleons in 
intermediate states to lie outside the Fermi sea. It should be stressed that 
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Figure 2.14: The diagram of the BBG expansion considered in Ref. |69] to 
obtain the high k and E behavior of the spectral function. 



Eq. ( I2.74P corresponds to the lp-2h diagram of Fig. 12. 141 which, through the 
G matrix, sums up all "ladder diagrams" where the BBG scattering matrix 
is replaced by the bare NN interaction. 

Eq. fl2.74p can be expressed in terms of the defect wave function, defined as 

1^12 >= 1^12 > -1012 > (2.76) 

where |\[^i2 > is the correlated two-body wave function, and |0i2 > the 
uncorrelated one. After a length algebra, one finds [53] 

P''^'{k,E) = ^^|^'^'"'c/A;cMfccMng^(fccM)n.e. (^^^^^ " ^^L/ + ^^o' 

(2.77) 

where 

P = (2.78) 
is the nuclear matter density and 

ko = ^Sm^v^ {e-EIH) . (2.79) 

It can be seen that the convolution formula fl2.64p is recovered. This repre- 
sents a robust validation of the convolution formula, as shown in Fig. 12.151 



57 



Chapter 2. The spectral function and the nucleon momentum 
distributions: two- and three-nucleon correlations 



strength functions 



5 



k - 



3 fm 



-1 



> 








100 



200 



300 



400 



E ( MeV ) 



Figure 2.15: The spectral function obtained from BBG theory plotted at 
/c = 3 fm~^ (diamonds). For comparison also shown are the spectral func- 
tion calculated in BBG theory with the free single particle spectrum (oc- 
tagons) and the one calculated by the convolution formula fl2.64p with the 
bare nucleon mass (full line), and the the pertinent effective mass (dashed 
line). After Ref. [69] • 

2.6 Many-body validation of the factorization 
of the nuclear wave function at high mo- 
menta 

In this section, we will show that the convolution formula results from a 
rather general property of the nuclear wave function. Let us recall that the 
convolution formula has been obtained by assuming that the nuclear wave 
function factorizes into the wave function of a correlated pair and the 
wave function of a core of the {A — 2)-nucleon system. In Ref. [70], the s- 
wave three-body Faddeev wave function obtained with MT/V [71] potential 
has been used to calculate the following quantity 



r = 



(2.80) 
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where 
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Figure 2.16: Eq. fl2.80p plotted versus the relative coordinate (12.811) . for 6 = 
90° and fixed values of y and y'. ( Upper panel) Green triangles: y = 3.0 fm, 
y' = 4.5 fm; black squares: y = 2.0 fm, y' = 5.0 fm; red circles: y = 5.0 fm, 
y' = 9.0 fm. {Lower panel) Red circles: y = 0.0 fm, y' = 1.0 fm; black 
squares: y = 0.5 fm, y' = 1.5 fm [70] . 
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|x| = |ri-r2| (2.81) 
is the two- nucleon relative coordinate, 



ri + r2 
^3 t: — 



(2.82) 

the CM coordinate, and 6 the angle between them. If the ratio fl2.80p is plot- 
ted versus |x| for fixed values of |y|, |y'| and 6, exhibits a constant behavior; 
this is evidence of the factorization of "^aIx, y, 0) in the variable x and y. As 
a matter of fact, it can be seen that writing 

^A{x,y,e) = f{x,e)h{y) (2.83) 

one gets 

^^{x^y^e) _ fix,e)h{y) _ %l-eo..t (2 84) 

^A{x,y',e) f{x,e)h{y') h{y') " ^'-'^^ 

The results presented in the upper panel of Fig. I2.16l clearly show that, when 
|y| is large and |x| small, indeed factorizes. On the contrary, as shown 
in the lower panel of Fig. \2.16\ for small values of |y| and small values of 
|x|, the nuclear wave function does not factorize. Thus the factorization of 
the nuclear wave function at high momenta is validated by many-body wave 
functions. 



2.7 Three-nucleon correlations 

2.7.1 Naive three nucleon correlation models 

We have seen that the high momentum and high removal energy components 
of the nuclear spectral function are due to few nucleon correlations (FNC), 
in particular in the region around 

FNC cannot however explain the regions at high values of k and large values 
of the excitation energy 

and the region at high values of k and low values of 
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■ Realistic P(k,E) (AV18)= 
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Figure 2.17: The realistic (black squares) spectral function of ^He [66j com- 
pared with the theoretical calculation performed within the convolution for- 
mula ^aMt) . at k = 3.5 fm-\ After Ref. 07]. 



This can be clearly seen in Fig. 12.171 In these regions, one has to consider the 
effects of three-nucleon correlations (3NC), i.e. configurations characterized 
by three nucleons which have comparable momenta, and share almost the 
full momentum of the nucleus 

ki + k2 + kg + Ka-3 = (2.88) 

i.e. 

ki + ka + kg ^ . (2.89) 

The naive 3NC model will be defined as follows: referring to Fig. 12. 18^ we 
will consider that the high momentum ki = k of nucleon '1' is completely 
balanced by the high momenta k2 and kg of particle '2' and '3', respectively, 
with the residual {A — 3)-nucleon system at rest. 

Within such a naive model, the residual system {A — 2) has momentum 

Ka-2 = kg = -(ki + ks) = -kcM (2.90) 
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Figure 2.18: The naive 3NC model: the high momentum ki = k of nucleon 
'1' is entirely balanced by the momentum k2 of nucleon '2' and ks of nucleon 
'3', with the system {A — 3) being at rest, i.e. with momentum K^_3 ~ 0. 
Thus, the momentum K.a-2 is given by the momentum ks nucleon '3'. 



and excitation energy 

2/i3,(A-3) A-2 2mN ^ ' 

given by the relative motion of particle '3' and the residual [A — 3)-nucleon 
system, with reduced mass 

A — 3 

/i3,(A-3) = . (2.92) 

The excitation energy E\_^ of the residual system {A — 1) can be written as 
follows 

E\ , = Mzh£ + (2.93) 

2/i2,(A-2) 2/i3,(A--3) 

i.e. by considering both the contribution arising from 2NC, given by Eq. 
fl2.55p . and the one due to 3NC, given by Eq. f l2.9ip . It can be clearly seen 
from Fig. 12. 18^ and it will be explained in more detail in §4.3[ that the 
assumptions of the naive model for three correlated nucleons refer only to 
that values of E\_^ given by Eq. (12.871) . 

The 3NC contribution to the spectral function P^{k, E), given by Eq. (I2.10p . 
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turns out to be [72] 

^ ^'^^ 2m;v(A - 1) V - 1) ; 2m;v(A - 2) ^"^^ 

X (5 (k + k2 + K^_2) (2.94) 



where A/3 is necessary to satisfy the normalization condition (12. 7p . It should 

soft 



be pointed out that Eq. f l2.94p is nothing but Eq. fl2.62p with ri^/^ replaced 



by n^c'lt- 



k, ^ 

' 2 




k3^-^ WK,_3-0 



Figure 2.19: First approximation to 3NC: the high momentum ki = k of 
nucleon '1' is entirely balanced by the momentum k2 = ka = — k/2 of nucleon 
'2' and '3', with the system {A — 3) being at rest, i.e. with momentum 
Ka-3 ^ 0. 



A first approximation to Eq. (I2.94p could be obtained by considering 

k2 = ks = -| (2.95) 

with 

Ka-2 = "2 = -kcM (2.96) 

and 

Ka-3 = (2.97) 

as shown in Fig. 12.191 

In this case, the "hard" CM nucleon momentum distribution will read as 
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follows 

nl^L\kcM) = S (kcM - (2.98) 
and the removal energy will reduce to 

E = eSI + ir-T7-^- (2.99) 



A-1 Am%j 



Integrating Eq. (12 .941) over the momentum k2 = — k — Ka-2, yields 



(3) _ A-3 Jl_ 
A-1 ArriN 



(2.100) 
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Chapter 3 

Inclusive electron scattering off 
nuclei at high momentum 
transfer and final state 
interaction effects: results of 
calculations 



Introduction 



Inclusive electron scattering off nuclei at high momentum transfer can pro- 
vide non trivial information on nuclear wave function. In particular, the 
kinematic region corresponding to Xbj > 1 is strongly affected by high mo- 
mentum and high removal energy components of the nuclear wave function 
arising from NN SRC. 

The Feynman diagram depicted in Fig. 13. II describes the inclusive A{e,e')X 
cross section which, as is well known, in one photon exchange has the follow- 
ing form 



dQo dv 



(3.i: 



where 



O? COS^ I 

4e^ sin^ | 



(3.2) 



is the Mott cross section, which describes the scattering from a pointlike 
nucleon, 



-9 



q2_zy2 ^4^^ sin2(^/2) >0 



(3.3) 
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Figure 3.1: Inclusive electron scattering off a nucleus A, in the one photon 
exchange approximation. In inclusive processes, denoted as A[e, e')X, only 
the scattered electron e' is detected in the final state; k, k', Pa and Px are 
the four-momenta of the incoming electron e, the scattered electron e', the 
target nucleus A, and the undetected particles, respectively; q is the four- 
momentum transfer such that — = = ~ where q = ki — k2 and 
u = ei — €2 are the three-momentum and the energy transfers, respectively. 



is the squared four-momentum transfer, 6 is the scattering angle and Wi^2 
are the nuclear structure functions. 

In plane wave impulse approximation (PWIA), shown in Figs, 3.2(a) and 



3.2(b), the latter are given by 
W,^{Q\u) = Y, f f dEPN{k,E) [CiW,''{Q^,u') + DiW^''{Q^,u' 

N=l 

(3.4) 

where % = {1,2}, are the nucleon structure functions, Cj and Di are 
kinematical factors whose explicit expression is given e.g. in [71], and, finally, 

= ^ (3.5) 

rriN 

where p is the four-momentum of the struck off-shell nucleon. 
According to the value of the invariant mass W produced by the interaction 
of the virtual photon with a nucleon in the nucleus, the inclusive process 
A{e, e')X is governed by the following two mechanisms: i) the quasi elastic 



66 



(a) (b) 




(c) (d) 

Figure 3.2: Processes contributing to the A{e, e')X cross section: (a) one- 
nucleon emission within the lA; (b) virtual photon absorption by a correlated 
NN pair within the lA; (c) single nucleon rescattering of a nucleon knocked 
out from shell model states; (d) elastic two-nucleon rescattering between the 
emitted nucleons of a correlated pair [73j . 



(qe) process, for which W = uin, and ii) the deep inelastic scattering (DIS), 
corresponding to W > ttin- The basic nuclear quantity that governs qe and 
DIS processes at xbj > ^ is the nucleon spectral function. 
The PWIA is only the first approximation to the calculation of the cross 
section, since the final state interaction (FSI) of the struck nucleon with the 
residual system {A — 1) may play an important role. As increases, the 
effects due to the FSI in the qe process are expected to decrease, whereas the 
DIS contribution increases. 

In this Thesis, we will focus only on the qe process. Two types of FSI in the 
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qe process, to be considered in what follows, are shown in Fig. 3.2(c) and 



3.2(d) 



3.1 The quasi elastic cross section 

In PWIA, the qe inclusive cross section describing the knock-out by the 
incoming electron of a nucleon from the nucleus A, can be written as 
follows [12] 



^2 (g, u) ^ 



J2 [dE f dkP^{k,E)aeN{q,^,k,E)6{u + MA-EN-EA-i] 



N=l 



(3.6) 



The argument of the energy conserving 5-function is 



u + Ma = sjm% + (k + q)2 + sjM*^_^ + k^ (3.7) 

where 

Ml_i = Ma-1 + E\_^ (3.8) 

is the mass of the excited residual system, and momentum conservation reads 
as follows 

q = p + P^_i = p-k. (3.9) 

Here z/ = ei — 62 and q = ki — k2 are the energy and three-momentum 
transfers, aeN is the elastic electron cross section off a moving off-shell nu- 
cleon with momentum = |k| and removal energy E, p and 'Pa-i are the 
momenta of the undetected struck nucleon and the final {A — 1) system; 
eventually, P^{k, E) is the nucleon spectral function discussed in the pre- 
vous chapter. From now on, for ease of presentation, only isoscalar nuclei, 
i.e. with P^{k,E) = P^{k,E) = P^{k,E), will be considered. 
After integrating over cos a = (k ■ q)/(A;g), Eq. (13.61) becomes 

a^{q,u) = 27r^ / dE kdkP^{k,E) 



{q,u,E) 

X aeAr(g, z^, k, E) 



du ^ 



kd cos a 



(3.10) 
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where the hmits of integration [75 



E 



\E 



A~l\ 



Ma^i +mN- Ma 



M*a-Ma 
(u + Ma) 



(3.11) 
(3.12) 



kcM -[q~ii^ + Ma)] J {Ma^i + E*^_,Y + K 



,2 

""CM 



MX 



(3.13) 



+ Ma) jfccA/ + [q-{i^ + Ma)] ^(M^_i + + A^^^^.j 



MX 



are imposed by energy conservation, where 



m: 



^^y + MAf-q" 



(3.14) 



(3.15) 



is the invariant mass, 




4(M^_i + E*_i)' 



N 



2M*, 



and tlie phase space factor 

du 



(3.16) 



(3.17) 



kd cos a ^/m% + q"^ + k"^ + 2kq cos a 

results from the dependence upon cos a of the nucleon energy in the final 
state in the 5-function. Eventually, 



O'eN 



O'Af 

El E2 



2\ 2 



(Ei + E^)' 



+ 



4 



(Fiat + E2n) 



e sin^ « (F^^ + rFl^) + -| (Fijv + Fajv) 



(3.18) 



is the electron-nucleon cross section for a relativistically moving nucleon, 
averaged over the polar angle, where Ei = qfi = q'^ — i^'^ , qf, = 

q2 _ - E2)\ and f = g2/(4 m^). 
At high momentum transfer, the quantity 



9z/ 



cos a 



(3.19) 
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depends very weakly upon k, so that it can be taken out of the integral and 
evaluated, e.g., at k = kmin- Therefore Eq. f l3.10p can be written in the 
factorized form 



Zaep + N 



du 


-1 


kd cos a 





F^{q,u) (3.20) 



where 



F^(g, z/) = 27r / dE kdkP^{k,E) (3.21) 

(q,v,E) 

is the nuclear structure function. 



3.2 The scaling function 

At high values of the momentum transfer, the rapid falloff of P^{k, E) with 
k and E allows the replacement Emax = kmax = +oo, and Eq. fl3.2ip can be 
written as 

r+oo r+oo 

F^{q,u) = 27i dE kdkP^{k,E) (3.22) 

(q,u,E) 

Therefore, the structure function depends upon q and v only through kmin, 

which is determined from the energy conservation (13. 7p . 

Let us replace the energy transfer v with a generic scaling variable 



Y = Y{q,v) 



(3.23) 



which is only required to be a function of q and v (and any arbitrary con- 
stant) so that, no matter with the specific form of Y , the cross section and 
the structure function can be expressed not in terms of the two canonical 
independent variables q and z/, but, without loss of generality, in terms of q 
and Y = Y{q, u). Correspondingly, a scaling function F^{q, Y) is introduced, 
which is nothing but Eq. ( I3.22p with u replaced everywhere by Y; if, under 
certain conditions, F^{q,Y) reproduces the asymptotic scaling function 



F^{Y) = 2n 



dE 



+ 00 



kdk P^{k,E) 



(3.24) 



Y-scaling is said to occur and, depending on the physical meaning of Y and 
F^(Y), various information on nucleons in nuclei could be obtained. 



70 



3.2. The scaling function 



Using Eq. (12.81) . the scahng function becomes 

F^(g,F) = 27r^ / naik) + 2TT dE kdkP^{k,E) 

(3.25) 

and it can be trivially cast in the form 

F\q,Y) = f\Y) - B\q,Y) (3.26) 
where the longitudinal momentum distribution 

r+oo 

f^{Y) =27r / kdkn^ik) (3.27) 

J\Y\ 

is integrated over all excited states of (A — 1), whereas the binding correction 

B^{q,Y) =2Tr / dE kdkP^{k,E) (3.28) 

on the contrary, is governed through kmin{q,Y, E) by the continuum energy 
spectrum of the final (A — l) system. The contribution arising from the latter 
strongly depends by the difference between Y and fcmm and, therefore, upon 
the definition of the former. In the Deuteron case, in fact 

E = E^in = 2.22 MeV (3.29) 

kmin{q,y,Emin) = \y\ (3.30) 

and thus 

i?^(g,r)=0 (3.31) 

F^(g,r) = r(F). (3.32) 

Therefore, when the binding correction can be neglected, as in the Deuteron, 
the quantities f^(Y) and n^{k) are linked by the relation 

1 df^(Y) 

SO that, if f^(Y) could be extracted from the experimental data, n"^(/c) could 
be determined. 

Unfortunately, in general, such an extraction is hindered by the presence of 

5^(g,r)^0 (3.34) 
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which leads to 

F^iq,Y)^nY) (3.35) 

and 

'df^(Y) dB\Y)' 



27tY 



dY dY 



(3.36) 



The factorization given in Eq. 03.26 p . in the asymptotic limit reads as 

pA^Y) = f^{Y) - B^{Y) (3.37) 

but again, unfortunately, owing to the presence of B^(Y), F'^iY) is not 
related to a momentum distribution, so that, in principle, the experimental 
longitudinal momentum distribution f^{y) and, consequently, n^^{k), cannot 
be extracted from the data. 

In both Eqs. 03.241) and 03.261) . the contribution arising from the binding 
correction depends upon the difference between Y and kmin{(l,Y, E), which 
could be minimized only by a proper choice of the scaling variable Y, such 
that kminil, Y, E) ~ \Y\. The resulting cross section 03.61) would depend only 
upon the nucleon momentum distributions, obtaining, by this way, a direct 
access to high momentum components generated by SRC [TB]. It is clear 
that the outlined picture can in principle be modified by the effects of the 
FSl of the knocked nucleon with the residual {A — l)-nucleon system. 



3.3 The final state interaction 



Owing to the decomposition rule 02. 8p of the spectral function P'^^k, E), the 
inclusive cross section can be written in the following form 



CTr, 



(3.38) 



where describes the transition to the ground and one-hole states of the 
(y4 — l)-nucleon system, and the transition to more complex highly excited 



configurations. The diagrams in Figs. 3.2(a) and 3.2(b) refer, respectively, to 



the contribution of and in PWIA, which is based upon the assumption 
that the reaction is well described by the exchange of a single virtual photon 
with a single nucleon, which does not interact with the remaining {A — 1)- 



nucleons; in particular, the diagram 3.2(a) represents the contribution arising 



from shell model configurations, whereas the diagram 3.2(b) mainly describes 
7* absorption by 2p-2/i configurations generated in the target ground state 
by NN correlations. It is well known that the PWIA sizably underestimates 
the cross section, both in light and heavy nuclei, at low values of z/, i.e. at 
XBj > 1, owing to FSl processes which arise from the interaction between 
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the knocked out nucleoli and the residual system. 

In this Thesis, following the approach proposed in Ref. [73], we take into 
account two different FSI processes, i.e.: i) the two-nucleon rescattering in 
the final state when the struck nucleon is a partner of a correlated pair, shown 



in Fig. 3.2(d), and ii) the single-nucleon rescattering of the struck nucleon 



with the mean optical potential of the residual {A ~ 1) system, depicted in 



Fig. 3.2(c) Let us discuss in more detail these two mechanisms. 



3.3.1 Two-nucleon rescattering 

The basic assumption underlying the convolution formula f l2.64p . is that two 
nucleons are locally strongly correlated at short separations, with their CM 
being apart from the spectator {A — 2)-nucleon system. 2NC in a nucleus are 
reminiscent of correlations in the Deuteron; indeed, as discussed in §2.2^ nu- 
cleon momentum distributions of light and complex nuclei, at A; > 1.5 fm~^, 
turn out to be the properly rescaled version of the Deuteron momentum dis- 
tribution. Therefore, in the 2NC region, the absorption of 7* by a correlated 
pair is expected to resemble the one in the Deuteron; if so, such a Deuteron- 
like picture of the initial state should be extended also to the final state by 



allowing the two nucleons to elastically rescatter, as depicted in Fig. 3.2(d) 
An important difference from the Deuteron case is that a correlated pair in 
a nucleus is bound and moves in the field created by the other nucleons. 



Within the above picture, the contribution of diagrams 3.2(a) 3.2(d) reads 
as follows 



= J2 f n^ik) [Z^aep + N^a^n] 6 {u + - Ek+,) (3.39) 

a<F •' 

= AaMott Yl J dkcMn^'M'ikcM)L''''W^j''' (3.40) 



NiN2=n,p ■ 



where 



Ep = \ ml + p^ (3.41) 



k'^ = MA-ViMA+\ea\-mMy + k^ (3.42) 

l^fj-f represents the (reduced) leptonic tensor, and W'^'^ is the hadronic tensor 
of a correlated pair, which can be written as follows 

/12 /3l2 /3l2 

X 5 (^u + k'cM - \/ {Mt ? + (kcAf + q)^) (3.43) 
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Figure 3.3: Inclusive cross sections at ~ 2 {GeV/cY versus the energy- 
transfer u. Calculations have been performed using the free nucleon form 
factors of Ref [77], the ccl presciption of Ref [75] for and the RSC 



potential [79] for the NN interaction. Dotted line: lA (Figs,3,2(a) ) + 3.2(b) 



dashed line: IA+ two-nucleon rescattering (Figs. 3.2(a) +3. 2(d) );dot-dashed 
line: contribution from inelastic channels estimated as in Ref. [80]. After Ref. 
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where is the nucleon current, fc^^ = Ma — \J M\_2 + kc-M; 1/^12 > is the 
relative wave function of the initial state of a correlated pair, and |/i2 > its 
continuum final state. Eq. fl3.40p is based upon the following assumptions 
on final and initial A-nucleon state: 

l^ii> ~ i{|/i2> iPcAf > |^L2>} (3-44) 
|^^> ~ ^{|/3i2>|xf2''>K-2>} (3.45) 

where ^ is a proper antisymmetrization operator, \xi'2^ > is the CM wave 
function of the initial state of a correlated pair, and |Pca/ > its plane wave 
final state. It can be seen from Eq. f l3.43p that medium effects on the 
hadronic tensor of the pair are generated by the energy conserving ^-function, 
in which the intrinsic energy available to the pair is fixed by its CM four- 
momentum, and, therefore, by the momentum distribution n^)^^ appearing 
in Eq. fl3.40p . Even if the CM motion is neglected by placing 

n^],f = 5(kcM) (3.46) 

medium effects still would remain through the quantity 

{k%j^)max = Ma- Ma-2 (3.47) 

which is related to the two-nucleon break up threshold, i.e. the binding of 
the pair. 

The inclusive cross section have been calculated in [75] for the Deuteron 
using the RSC NN potential [TS], taking into account the rescattering in 
S, P and D partial waves; then, using the same two-nucleon amplitudes 
< PulJf^^ + j/^^|/i2 >, the cross section aj^ have been computed for complex 
nuclei. The results are shown by the dashed lines in Fig. 13.31 it can be 
seen that, at 1.3 < xbj < 2, the process of two-nucleon rescattering brings 
theoretical predictions in good agreement with experimental data. The most 
striking aspect of these results is that the same mechanism which explains 
the Deuteron data, does the same in a complex nucleus, provided the A 
dependence due to n^M^^ and fc^-M (clearly exhibited in Fig. 13. 3p is properly 
considered. It should be pointed out that these results hold for the whole set 
of kinematics considered in Refs. [131 ED IS21 1^ . 

3.3.2 Single nucleon rescattering 

The two-nucleon rescattering is not able to describe the experimental data 
at x^j ^ 2. This fact is not surprising because, at xbj ^2, more than two 
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nucleons should be involved in the scattering process. This process can be 
simulated by considering the motion of the nucleon, knocked out from shell 
model states, in the optical potential generated by the ground state of the 
(A— l)-nucleon system. Within such an approach, corresponding to diagrams 
3:2(i0|3.2(c)| Eq. (Km> becomes 



^ Vopt 



(3.48) 

resulting from the eikonal approximation for the nucleon propagator [8^ . 
The optical potential Vopt can be cast in the following on-shell form 

Vopt = -p Vn ctnn ^ (3.49) 

where p is the nuclear density, vn is the nucleon velocity, a^N is the total 
NN cross section and a^N is the ratio of the real to the imaginary part of the 
forward NN scattering amplitude. The imaginary component of Eq. f l3.49p 
describes inelastic processes leading to excitations of the residual system, 
which, in the high energy regime, mainly correspond to secondary nucleon 
emissions. 

It has be pointed out [73], that the two- nucleon rescattering is not included 



in the process described by diagram 3.2(c) As a matter of fact, indeed, the 



two-nucleon rescattering is not a multiple scattering process, i.e. it does not 
contribute to an optical potential, and, particularly, is independent of nuclear 
density and does not produce any absorption of the outgoing flux. 
However, treating FSl at x^j > 1 in terms of on-shell potentials is not justi- 
fied [HI Ul ESI [86]. Indeed, the struck nucleon, with momentum 

p'2~(z/ + m^-E)2-(k + q)2 (3.50) 

can be either on-mass-shell, i.e. with p'"^ = or off-mass-shell, namely 
with p'"^ ^ m^, depending on the values of k and E. Initial configurations 
with k < kmin give rise to an intermediate off-mass-shell virtual nucleon, so 
that the use of an on-shell optical potential is unjustified, and can hardly be 
reconciled with the fact that rescattering amplitudes are expected to decrease 
with virtuality, because an off-shell nucleon has to rescatter within short 
times. Therefore, in order to take into account off-shell effects, we have to 
include in Vopt a suppression factor of the type 

Vopt = -^pVN (^NN (« + ajviv) e^'^l*^'~p"l . (3.51) 

The differences between the use of an on-shell or off-shell potential are shown 
in Fig. 13.41 The parameter 6 appearing in Eq. fl3.5ip is the same for all 
kinematics considered. 



76 



3.4. Inclusive cross sections: results of calculations 




0.6 0,8 1 1-2 1.4 

V (GeV) 



Figure 3.4: The cross section (Tq versus the energy transfer z/. Dotted line: Eq. 



(Km (diagram [MJ^; dashed line: Eq. (K^ (diagrams [3:2(iO] + |3.2(c)D 
calculated using the on-shell potential fl3.49p : solid line:Eq. fl3.48p (diagrams 



3.2(a) + 3. 2(c) ) calculated using the off-shell potential f l3.5ip . After Ref. [73] . 



3.4 Inclusive cross sections: results of calcu- 
lations 

The results of our calculations, including the contributions from all diagrams 
depicted in Fig. 13. 2[ are shown in Figs. 13.51 and 13.61 It can be seen that 
the PWIA overestimates the experimental data at xbj > 1-5, whereas the 
inclusion of the FSI produces a good agreement between theoretical calcula- 
tions and experimental data. In the region 1.5 < xbj < 2, the FSI is mainly 
due to the two-nucleon rescattering whereas, at xbj > 2, the contribution 
from the optical potential, which mocks up three-nucleon correlation effects, 
becomes important. Other results are shown in Figs. 13.71 and l3.8[ leading to 
the same conclusions. Results of the same quality were previously obtained 
in Ref. [HI HSl [73] . 
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Figure 3.5: The experimental inclusive cross section ^^Fe{e, e')X ^83j vs. the 
energy transfer z/, compared with theoretical calculations which include SRC 
and FSl. Dot-dashed line: PWIA, Eq. (K^ : dashed line: PWIA + FSl of 
the correlated struck nucleon with the correlated partner, Eqs. f l3.39p and 
f l3.40p : solid hne: the same as dashed red line plus the FSl of the shell model 
struck nucleon with the mean field of the residual {A — l)-nucleon system, 
Eqs. fIXlg]) and f MI]) . After Ref [SP. 
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Xg= 1.94 1.6 1.42 1.19 1.07 




1.00 1.25 1.50 1.75 2.00 



v[GeV] 

Figure 3.6: The same as in Fig. 13. 5[ for different fcinematic conditions. After 
Ref. 
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Figure 3.7: The same as in Fig. 13.51 Experimental data from Ref. [Ti] . 



3.4. Inclusive cross sections: results of calculations 



E=5.766 GeV 
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Figure 3.8: The experimental inclusive cross section ^^Fe{e, e')X [S7| vs. the 
energy transfer p, compared with our preliminary calculations which include 
SRC and FSl. Dot-dashed line: PWIA, Eq. fl3^ : solid line: PWIA plus 
the FSl of the shell model struck nucleon with the mean field of the residual 
{A — l)-nucleon system and the FSl of the correlated struck nucleon with the 
correlated partner, Eqs. fl3.48l) and fl3.4UI) . 
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Chapter 4 

A novel approach to scaling 
phenomena in inclusive 
scattering: mean field, 
correlations and proper scaling 
variables. 



Introduction 

The inclusive cross section can be analyzed in terms of scaling function and 
scaling variable. 

The scaling variable Y arises from the energy conservation law which, for an 
inclusive process, reads as follows 

u + Ma = \J {Ma-1 + El_^f + + ^Jml + (k + q)2 (4.1) 

being 

Ma = mass of the target nucleus (4.2) 

Ma-i = mass of the residual nucleus {A — 1) (4.3) 

E*A_^ = intrinsic excitation energy of the residual system (4.4) 

niN = mass of the knocked out nucleon (4.5) 

As mentioned in Chapter [21 a proper choice of the scaling variable Y could 
minimize the contribution arising from the binding correction, allowing a 
direct link between the scaling function and the nucleon momentum distri- 
bution. Within such an approach, the new inclusive cross section for electron 
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scattering off nuclei can be written only in terms of nucleon momentum dis- 
tributions, so that inclusive scattering becomes a powerful tool to investigate 
NN SRC in nuclei. In what follows, it will be illustrated in detail how the 
dependence of kmin upon E\_^ gives rise to the binding correction, and how 
to a different definition of the scaling variable Y. 



4.1 The mean field scaling variable 

The final state interaction of the struck nucleon invalidates the PWIA but, 
in spite of that, an approach was developed in the past to reduce the effects 
from both the binding corrections and FSI [HI [75]; the traditional approach 
to Y-scaling is based upon the traditional scaling variable Y = y, obtained 
by placing 

k = \y\ (4.6) 

cos a = — r— ^ = 1 (4.7) 
kq 

E\_, = (4.8) 
in the energy conservation law given by Eq. (14. obtaining 

u + MA = ^Ml_,+y^ + ^m% + {y + qy (4.9) 
whose solution is [12] 



-qA ± -^^2^2 _ ^ ^^^2 _ ^2 

with 



(4.10) 



Ml ={u + MAy-q' (4.11) 
and 

A = + Ml_i - m% . (4.12) 

Within such an approach, y represents the minimum longitudinal momentum 
of a nucleon having the minimum value of the removal energy 

E = Ermn + = ^min = TRn + Ma-1 - Ma ■ (4. 13) 

At high values of g, one has 

lim krmn{q,y,E) = kZniV , E) = \y-{E- Ermn)\ (4.14) 
g— >-oo 
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so that, when E = Emin, 

kZn{y,E) = \y\ (4.15) 
B^{q,y) = 27r / dE kdkPf{k,E) = (4.16) 

^ ^min ^ \ y\ 

and thus the scaling function reduces to 

f+oo 

F\y) = f\y)= kdkn\k) (4.17) 
J\y\ 

explicitly showing scaling in y. 




-600 -400 -200 

y [MeV/c] 



Figure 4.1: The ratio of the binding correction B^{q,y) (Eq. fl3.28p ) to the 
scaling function F^(g, y) (Eq. fl3.2ip ) for ^He (open dots) and (full dots), 
calculated using the scaling variable y. After Ref. [76] . 

Unfortunately, this occurs only in the Deuteron, whereas in the general case. 
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for A > 2, the excitation energy E\_i of the residual system is different from 
zero and 

E = Emin + E\_^ > Emin (4-18) 

leading to 

B^{q,y)>Q (4.19) 

and thus to the relation F^{y) ^ f^{y), given by Eq. fl3.26p . 
To illustrate the relevant role played by the binding correction in the tradi- 
tional approach to Y-scaling, the ratio 

F^{q,y) ny)-B^{q,y) ^ " ^ 

which represents the deviation of the scaling function F^{q,y) from the lon- 
gitudinal momentum distribution f^{y), is shown in Fig. 14. ![ plotted versus 
the scaling variable y. It can be seen that, at high values of \y\, the effects 
from binding are very large whereas, at low values of \y\, binding effects can 
be neglected. 

Moreover, the experimental scaling function 

FfA<i,y) = i (4.21) 



plotted versus the scaling variable y, as shown in Fig. 14. 2[ confirms that 
the scaling function strongly differs from the longitudinal momentum distri- 
bution, and therefore does not exhibit any proportionality to the Deuteron 
scaling function f^{y). Therefore it should be pointed out that, when ex- 
pressed in terms of y, a comparison between experimental and theoretical 
scaling functions requires the knowledge of the nucleon spectral function, 
generated by the main role played by B'^{q,y). 

Moreover, the experimental scaling function exhibits a strong q dependence 
owing to the FSI and binding effects, and differs from the asymptotic scaling 
function F^{y). The latter, however, has been obtained in Ref. [HI |75] by 
extrapolating to g — ?■ oo the available values of F^{q,y), on the basis that 
FSI can be represented as a power series in 1/q, and dies out at large g^, a 
conclusion that has been reached by various authors (see e.g. [5U]). 
It is therefore the dependence of kmin upon E\_^ that gives rise to the bind- 
ing effect. This is an unavoidable defect of the usual approach to Y-scaling, 
based on the scaling variable y; except for the trivial case of the Deuteron, 
in fact, in a complex nucleus the final {A — l)-nucleon system can be left in 
all possible excited states, including the continuum but, by definition, the 
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Figure 4.2: The experimental scaling function F^p{q,y) of ^He, ^^C, and 
^^Fe obtained from the experimental data of Refs. P^IBS]. The longitudinal 
momentum distributions (Eq. f l3.27p ) of (full line), ^He (long-dashed), 
^^C (dashed) and ^^Fe (dotted) are also shown. After Ref. [76] . 



traditional scaling variable y can only be identified with the longitudinal 
momentum of weakly bound, shell nucleons {E\_^ ~ — 20MeV). The 
longitudinal momentum for such nucleons is very different from the strongly 
bound, correlated nucleons {E\_^ ~ 50 — 200 MeV), and this explains why, 
at large values of \y\, the scaling function is not related to the longitudinal 
momentum of strongly bound correlated nucleons, whose contributions al- 
most entirely exhaust the behavior of the scaling function. 
As stressed in Refs. [TBI IIZl US], to establish a global link between exper- 
imental data and longitudinal momentum components, one has to conceive 
a scaling variable that could equally well represent longitudinal momenta of 
both weakly bound and strongly bound nucleons, so that the binding correction 
could be minimized. 

The experimental longitudinal momentum distribution f^{y) has thereby 
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been obtained by adding to F^{y) the binding correction B^{y) evaluated 
theoretically, as shown in Fig. I4.3[ and n^^{k) has been obtained by Eq. 
fl3.36p . Such a procedure affects the final results in terms of large errors on the 
extracted momentum distributions, particularly at large values of k] in spite 
of these errors, the extracted momentum distributions, at k > 1.5 — 2 fm~^, 
turned out to be larger by orders of magnitude from the prediction of mean 
field approaches, and in qualitative agreement with realistic many-body cal- 
culations that include SRC. 

In order to make the extraction of f^{y) as independent as possible from 
theoretical binding corrections, in Ref. [TB] is thus necessary to introduce an- 
other scaling variable, incorporating relevant physical dynamical effects left 
out in the definition of y. 




Figure 4.3: The longitudinal momentum distribution for (dotted line), 
"^He (full line) and ^^Fe (dashed line) corresponding to a parametrization 
obtained in Ref. [T7] . 
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4.2 Two-nucleon correlation scaling variable 



2NC are defined, as previously explained in §2.41 as those nucleon configura- 
tions shown in Fig. 14.41 1!3]: momentum conservation in the ground state of 
the target nucleus 

A 

(4.22) 

i=l 

is almost entirely exhausted by two correlated nucleons with high momenta, 
the {A — 2)-nucleon system acting mainly as a spectator, moving with very 
low momentum. The intrinsic excitation energy of the (yl — l)-nucleon system 

\< A \A > 

\ 




K^_l\ K^_2 = i 



I 



Figure 4.4: 2NC correlations in a nucleus A: the high momentum ki = k of 
nucleon " 1" is almost completely balanced by the momentum k2 ~ — k of the 
partner nucleon "2", whereas the residual system moves with low momentum 
Ka-2- Momentum conservation is kj = ki + k2 + K^_2 = 0. 

is in this case 
which becomes 

in the naive 2NC model, i.e. the model based upon the assumption Ka_2 = 0. 
Since high excitation states of the final {A — l)-nucleon system are generated 
by SRC in the ground state of the target nucleus, the traditional (mean 
field) scaling variable y does not incorporate, by definition, SRC effects, for 
it is obtained by placing E\_^ = in the energy conservation law (14. ip . 
Motivated by this observation, in Ref. [16], a new scaling variable, Y = 
Ucw = 2/2 has been introduced, by setting 

k = \y2\ (4.25) 
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cosa = ^V^ = l (4.26) 
kq 

E\,^ =< E\,^{k) >2NC (4.27) 
in Eq. (14.11) . which in this case reads as follows 



(4.28) 

The inclusion of the term 

< E\_,{k) >2^c= ^ j P^r,c{k.E\^,)E\_,dE\_, (4.29) 

makes the scaling variable ?/2 to properly depend upon the momentum depen- 
dence of the average excitation energy of (A — 1), generated by two-nucleon 
correlations. Here 

E\_, = E- Elll (4.30) 

where 

eIH = Ma-2 + 2mN - Ma (4.31) 
is the threshold energy for two-particle emission. 

The quantity in Eq. f l4.29p has been calculated using a realistic spectral 
function for nuclear matter and ^He. The results are presented in Fig. 14. 5[ 
where they are compared with the prediction of the spectral function of the 
few nucleoli correlation (FNC) model of Ref . [00] , according to which 

2 



A - I 2mN 



(4.32) 



where Kcm is the CM momentum of a correlated pair. In view of the 
very good agreement between the FNC model and the exact many-body 
results for nuclear matter and ^He, the former has been used to calculate 
< E\_^{k) >2NC for nuclei with 3 < A < 00. 
The values shown in Fig. 14.51 can be interpolated by 

<E\_,ik)>2NC=^^T^ + bA-CA\k\ (4.33) 

where 
and 



= \r^N + k^-mN (4.34) 



JdkcMtk^P,{k,E) 

= riv — Pfu — ^-35 

j dkcMPi{k,E) 
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J dkcM k ■ kcM Pi {k, E) 



result from the CM motion of the pair. 




Figure 4.5: The average value of E\_^{k) [Eq. fl4.29p ] calculated for nuclear 
matter with the spectral function of Ref. [63l El] (open dots), and for '^He 
with the spectral function from the Pisa wave functions [90] (full dots). The 
full lines are obtained with the spectral function of the few-nucleon correla- 
tion model of Ref. [nOj . 

The values of 6a and ca used in our calculations are listed in Table 14.11 
For a complex nucleus and not too large values of 1/25 a solution for Eq. f l4.28p 
can be written as _ 

y^--hw:.^^^^- (4.37) 

where 

pa = t^ + Md (4.38) 
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Nucleus 


bA (MeV) 


CA 




-2.94 


-0.03 


NM 


37.3 


0.04 



Table 4.1: Values of the parameters appearing in Eqs. fl4.35p and fl4.36p for 
complex nuclei and nuclear matter. 



Md = 2mjv - - bA+ < Egr > (4.39) 

q = q + CAiyA (4.40) 

Wl = -q^ = Ml + 2z/Md - Q' (4.41) 

In order to counterbalance the effects of < E^_^{k) >2nc at low \y2\, in the 
definition of 

Ml_i = Ma+ < E\_,{k) >2NC -<Egr> (4.42) 
has been added the value < Egr >, fixed by the Koltun sum rule (see [TB| 

In the Deuteron case 

y^ = y = -l + ^^^D-^rn% (4.43) 

with 

ujj = u + Md (4.44) 

and where 

= un^ - = Md^ + 2uMd - (4.45) 

is the Deuteron invariant mass. 
For small values of I1/2I, such that 

sjyl + m]^ - +hA- CA\y2\ « Egr > (4.46) 

the variable y, representing the longitudinal momentum of a weakly bound 
nucleon, is recovered. 

Therefore y2 effectively takes into account the k dependence of E\_^, both 
at low and high values of y2, and interpolates between the correlation and 
the single-particle regions; it can be interpreted as the minimum longitudi- 
nal momentum of a nucleon that, at high values of y2, has removal energy 
< E\_-^ >2NC o.iT'd is partner of a correlated two-nucleon pair with effective 



A-2 
A-1 
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mass Mu- 
lt should be pointed out that, in our calculations, the fourth-order equation 
resulting from Eq. (I4.28P has been solved exactly; this, together with the 
relativistic extension of the definition of the mean excitation energy, is nec- 
essary to extend 1/2 to high values. 




-600 -400 -200 

[MeV/c] 



Figure 4.6: The ratio of the binding correction B^{q, 1/2) (Eq. f l3.28p ) to the 
scaling function F^{q,y2) (Eq. fl3.2ip ) for ^He (open dots) and ^^C (fuU 
dots), calculated using the scaling variable 1/2 • After Ref. [76] . 

Within such an approach to Y-scaling, the effects due to the binding are 
strongly suppressed, as clearly illustrated in Fig. 14. 6[ where the ratio given 
by Eq. fl4.20p vanishes in the whole region of y2 considered; the main feature 
of 1/2, in fact, is that 

{q,y2,E) c^\y2\ (4.47) 

which leads to 

B^{q,y2)-0 (4.48) 
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with the following two relevant consequences: 

1. the relation 

r+oo 

F\q, y^) ^ f\y^) = k dk n\k) (4.49) 

J\y2\ 

holds, thus, plotting the data in terms of ?/2 can provide a direct access 
to the nucleon momentum distributions, and so to SRC in nuclei; 

2. as mentioned in § 2.3, many-body calculations show that at high mo- 
menta, k > 2 fm~^, all nucleon momentum distributions are simply 
rescaled version of the Deuteron one, i.e. 

n^(A;) ~ C^n^ik) (4.50) 

where is a constant; as a consequence, one would expect that, 
at high values of \y2\, F^{q,y2) will behave in the same way in the 
Deuteron and in complex nuclei, so that, accordingly, 

ny2)^C^f''iy2)- (4.51) 

At the same time, on the contrary, at low values of \y2\, F^{q,y2) 
should exhibit an A dependence generated by the different asymptotic 
behavior of the nuclear wave functions in configuration space. 

This is fully confirmed by Fig. 14.71 where the scaling function F^{q, 2/2) ob- 
tained from available experimental data on ^He, ^^C and ^^Fe, is plotted 
versus the scaling variable 1/2; it can be seen that, at high values of \y2\, 
F^{q,y2) scales exactly to f^{y2)- This is even better demonstrated in Fig. 
14.81 where the scaling function F^(g, 7/2) is plotted versus at fixed values 
of y2- By this way, the scahng behavior of F^{q,y2) is better illustrated. In 
the same Figure, in order to analyze more quantitatively the scaling behavior 
of F'^i^q, 7/2), the latter has been plotted together with the theoretical scaling 
function for A = 2, calculated in PWIA (solid line), and taking FSI into ac- 
count (dashed line). The left panel clearly illustrates that, due to FSI effects, 
scaling is violated and approached from the top, and not from the bottom, 
as predicted by the PWIA. However, scaling violation seems to exhibit a 
dependence which is very similar in Deuteron and in complex nuclei. This is 
illustrated in more details in the right panel of the Figure, which shows the 
quantity 

F^(Q^l/2)/C-^^F^(g^l/2) (4.52) 

where the A dependent constants have been chosen so as to make the 
experimental scaling function for a nucleus A to coincide as much as possible 
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[MeV/c] 



Figure 4.7: The experimental scaling function F^p{q,y) of '^He, ^^C, and 
^^Fe obtained from the experimental data of Refs. [HI [88]. The longitudinal 
momentum distributions (Eq. f l3.27p ) of (full line), ^He (long-dashed), 
^^C (dashed) and ^^Fe (dotted) are also shown. After Ref. [7B] . 

with the Deuteron scaling function F^i^Q"^ ,1/2). It clearly appears that the 
scaling functions of heavy and light nuclei scale to the Deuteron scaling 
function; moreover the values obtained for turn out to be in agreement, 
within the statistical errors, with the theoretical predictions of Ref. [6], as 
well as with the experimental results on the ratio [3] 



shown in Fig. 14. 9j it is also important to stress that, although FSl are very 
relevant, they appear to be similar in Deuteron and in a nucleus A, which is 
evidence that, in the SRC region, FSl are mainly restricted to the correlated 
pair. 
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Figure 4.8: The scaling function F^{Q'^,y2) from the lower panel of Fig. 
14.71 plotted vs. at fixed values of y2 {'^ He- asterisks , ^'^C -triangles, ^^Fe- 
squares). In the right panel the data for ^He, ^^C and ^^Fe have been 
divided by the constants = 2.7, C^^ = 4.0 and C^^ = 4.6, respectively. 
The theoretical curves refer to and represent the PWIA results [full) 
and the results that include the FSI (dashed), both obtained with the AV18 
interaction [20]. Scaling variables are in MeV/c. After Refs. jlTj 176] 
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Figure 4.9: Eq. f l4.53p vs. the scaling variable X = xbj for A = 56 at six 
different values of Q^. The average is given for each frame. To the right 
of the vertical dashed line those data correspond to a final state less 50 MeV 
greater than the Deuteron rest mass. The solid line is a calculation based on 
the nuclear spectral function of Ref . [21] . After Ref . |3] . 

4.3 Three- nucleon-correlation scaling variable 

Three-nucleon correlations also contribute, in principle, to the inclusive cross 
section for A > 3. 3NC correspond, as previously explained in §2.7^ to those 
three-nucleon configurations in which the high momentum ki = k of nucleon 
"1" is almost entirely balanced by the momenta k2 and ka of nucleons "2" 
and " 3" , respectively. 

The excitation energy of the [A — l)-nucleon system is given in this case by 
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(Cf. Eq. I^M)) M 

E* 



(k2 - kg)' ^ ^ - 3 [(k2 + ks) - 2K^„3] 



rriN 



A-1 



4m 



TV 



(4.54) 



which shows that, whereas 2NC are directly linked to high values of excitation 
energies 

{A - 2) P 



7-1* 



(4.55) 



{A - 1) 2mN 

high momentum components due to 3NC may lead both to low and to high 



values of E\_^, as shown in the examples of Figs. 4.10(a) and 4.10(b) 
respectively. 



k 



k. 



K,_3 - 



Si^k.-k 



(a) 



(b) 



Figure 4.10: Two types of 3NC configurations which are present in the spec- 
tral function of a nucleus A; they correspond to: (a) high momentum k and 
low removal energy E, and (b) high momentum k and high removal energy 
E. 



In the configuration of Fig. 4.10(a) the momentum ki = k of nucleon 
"1" is almost entirely balanced by nucleons "2" and "3", with momenta 
k2 ~ kg ~ — k/2, and one has 

A-3 e 



rp* 



A-1 AniN 



In the configuration of Fig. 4.10(b) 



with 



k2 = 



cos6' 



— cos , 
2 V2 



(k2 ■ ks) 

ik2ks) 



(4.56) 



(4.57) 



(4.58) 
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■ Realistic P(k,E) (AV18)= 




I I I I [ I ] 
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E [MeVl 

Figure 4.11: The spectral function of ^He vs. the removal energy at 
k = 3.5 fm~^ [66]. The realistic spectral function corresponding to the Pisa 
wave function (squares) is compared with the FNC model (full line) given by 
Eq. (IZMj) . at = 3.5 fm-\ After Ref. |17]. 



and E\_^ could be very large when k2 and are large. 

Let us investigate the presence and relevance of 3NC configurations in the 

spectral function of the 3N system, for which the Schrodinger equation has 



been solved exactly. When A = 3, 3NC of the type shown in Fig. 4.10(a) 
lead to E* = (cf. Eq. (gSSD). 

In Fig. 14.111 as already illustrated in §2.71 the realistic spectral function 
of ^He obtained [66j using the Pisa wave function [2D] corresponding to the 
AV18 interaction ^20j (full squares), is compared with the predictions of the 
FNC model (solid line) given by the convolution formula (12.641) [47J. It can 
be seen that the FNC spectral function reproduces the exact one in a wide 
range of removal energies (50 ^ E < 250 MeV), but fails at very low and 
very high values of E, where the effects from 3NC are expected to provide 
an appreciable contribution. 



can 



It is clear from Fig. 14.111 that 3NC of the type shown in Fig. 4.10(b) 
hardly be present at A; < 3.5/m~^ and E < 300 MeV", so that it is legitimate 
to ask ourselves whether these 3NC can show up in available experimental 
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data. To answer this question, let us now consider the maximum value of the 
removal energy achieved in the experiments, which is also the upper limit of 
integration in Eq. fl3.2ip . i.e. 

EmaM. ^) = VFTMa^^ . (4.59) 
In Fig. 14.121 we show the value of Emaxig,^) plotted versus the Bjorken 
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Figure 4.12: The maximum value of the removal energy Emax (Eq. fl4.59p ) 
available in inclusive q.e. scattering off ^He plotted vs. xbj, at increasing 
values of shown in the inset. After Ref. j47] . 



scaling variable in the region 1 < XBj < 3, in correspondence of a set of 



values of u and q typical of available experimental data on ^He 
seen from Figs. 14.81 and 14.111 that, in the region 2 < xbj < 3, 
configurations of the type shown in Fig 



It can be 
only 3NC 
can contribute to present 
A{e, e')X kinematics; for this reason we will consider, for the time being. 



4.10(a; 
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only this type of 3NC. 

Therefore the 3NC scaling variable F = ^3 is obtained by placing 

k = Ivsl (4.60) 
k • q 

cos a = — — = 1 (4.61) 
kq 

E\_^=<E\_^{k)>^NC (4.62) 
in the energy conservation (14. ip . which becomes 



u + Ma= ^J{Ma-3 + 2m^+ < E\_^{k) >3ivc)' + yl + \J m]^ + {Vz + q? 

(4.63) 

where the excitation energy of the residual system is (cf. (Eq J4.56l) ) 

<^-W>'--1^3^ («4) 



and corresponds to the 3NC configuration of type 4. 10 (a 



Even in this case, Eq. (14.631) has been solved exactly in our calculations. 

4.4 Domain of existence of the three scaling 
variables 

In the previous sections, we have obtained three different scaling variables, 
y2 and y^, by placing different values of E^_^ in Eq. (14. ip . namely E\_^ = 0, 
E*A^i = <E\ _i{k) >2NC, and =< E*^^^{k) >3nc, respectively. 
In Fig. I4.13[ the values of y, y2 and 7/3 are plotted versus xbj, in the case of 
^^Fe, for a fixed value of Q^. It should be pointed out that the magnitude of 

XB, = = A (4.65) 

cannot be larger than the number of nucleons in a given nucleus, as clearly 
results from the definition of the invariant mass 

_ _q2 j^m\^ 2Mav = (-1 + — ^ I +Ml>Ml (4.66) 

V XTUnJ 

which leads to 

XBi< — ^A. (4.67) 
rriN 
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Figure 4.13: The scahng variables 1/2 and vs. xbj for the nucleus ^^Fe. 
After Ref. [H]. 

This result is really important in order to explain the different behaviors 
of the three scaling variables which, as already discussed in the previous 
sections, differ only in the definition of E\_^ used in Eq. (14. ip . Let us 
analyze this point in more detail. 
The 2NC scaling variable fl4.37p . requires that 

Wl - Am\ > . (4.68) 
By approximating ~ 2m7v, and owing to Eq. fl4.65p . one gets 

- Q2 + 2 ^ > 1 (4.69) 

X 

and thus 

(l--\>l (4.70) 
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which means 



XBj < 2 . 



(4.71) 



The scaling variable y2 is therefore defined only for xbj < 2. Indeed ?/2 
represents 2NC in heavy nuclei resembling the ones acting in Deuteron and, 
by placing A = 2 in Eq. f l4.67p . the same domain of existence is recovered. In 
the same way, the domain of the other two scaling variables y and y^, turns 
out to exist in the range of x^j < A and xbj < 3, respectively. Indeed, the 
scaling variable y describes the mean field configuration, whereas the scaling 
variable y^ describes 3NC as in ^He; thus, by placing A and A = 3 in Eq. 
( I4.67p . the same domains are recovered. 

The different values of E'^_^ used in Eq. fl4.ip . impose therefore different 
limits in the domain of existence of the three scaling variables used in our 
calculations, due to the different configurations of nucleons in nuclei they 
refer to. 

4.5 A new approach to the treatment of the 
inclusive cross section 

Let us recall, as already discussed in §3.H that, at high Q^, the calculation 
of the quasi elastic inclusive cross section depends upon the scaling function 



whose calculation requires, due to the presence of the nucleon spectral func- 
tion P{k, E) (Eq. fl2.ip ). the knowledge of the entire energy spectrum of the 
{A — l)-nucleon system. Owing to the decomposition rule fl2.8p . the contri- 
butions from different final nuclear states could be explicitly separated out, 
writing 



F\q,y)=2ix 




kdk P^{k,E) 



(4.72) 



F^iq,y) 



foiy) + E2{y,q) + F3{y,q) 



(4.73) 



where 




(4.74) 



describes the shell model contribution 



F,^{q,y) = 2n 




kmin{y,q,E) 



k dk P2Ncik, E) 



(4.75) 
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section 



the 2NC pair contribution and, eventually, 

/•H-oo p+oo 

F,^iq,y) = 27T dE k dk P,r,cik, E) (4.76) 

the 3NC contribution. 

As already pointed out in the previous sections, the scaling variables y, y2 
and 2/3 effectively take into account the energy -E^„i of the residual system; 
by this way, the effects due to the binding correction f l3.28p are strongly 
suppressed, and a direct link between the scaling function and the nucleon 
momentum distributions can be established. 

The original idea of our novel approach [47j to A{e, e')X processes, is based 
upon the replacement of Eq. f l4.73p by the following one 

F;tJy,q)^fneJy) = h{y) + f2{y2) + h{yz) (4.77) 

with /o(y) given by Eq. fl4.74|) . and Eqs. ( I4.75P and ( I4.76p replaced by 



f^{y2) = 2n / n^ik) kdk (4.78) 
J\y2\ 



and 

"+00 



r+oo 

fsiys) = 2n / 4{k) kdk (4.79) 

respectively. It should be pointed out that y2 = 2/2(0', y) and 2/3 = y3{q,y), 
with y = yjq u). 

In Eq. 04.78 p . we use the nucleon momentum distribution 

^2 (^) = j dkcM iireii^ + kcM) n^°//(kcj\/) (4.80) 

which describes the virtual photon absorption by a 2N correlated pair and, 
in Eq. (14. 79 p . the nucleon momentum distribution 

n^{k) = j dkcM n.e/(k + kcM) n1,lf{kcM) (4.81) 

which corresponds to the contribution arising from 3NC; for shell model 
nucleons the usual nucleon momentum distribution (see Eq. f l2.42p ) 

^oik) = -^ E ^"^-(^) (4.82) 



AttA 
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defined in §2.31 is adopted. 

These three quantities should satisfy the normalization condition fl2.19p . i.e. 

/•oo i*co roc 

/ k^dkn^{k)+ k'^dkn^{k)+ k^dknf{k) = l. (4.83) 
Jo Jo Jo 

Due to the difiiculties in the calculation of the nucleon momentum distri- 
bution fl4.8ip . until now we have considered only the independent particle 
shell model and 2NC components in Eq. f l4.77p . We have thus compared 
the longitudinal momentum distribution fl3.27p . where E = E^in, with Eq. 
f l4.77p . where E = Emm < E\_^{k) >2nc^ and with the exact calculation of 
Eq. f l4.73p in terms of the spectral function, finding, as shown in Fig. I4.14[ 
a good agreement between the two approaches. 

It appears, therefore, that inclusive cross sections can be calculated only by 
using momentum distributions, provided the excitation energy of the {A — 1) 
system is effectively taken into account in the lower limit of integration in k, 
i.e. by using the scaling variable y2- 

The same conclusions can be reached also by writing the scaling function in 
a more general form, viz. 

F'^{q,Y)= dE d^kP^{k,E). (4.84) 

The scaling variables y2 and y-^ allow us to write 

P^{k, E) = n^{k) 6{E - E„,„) (4.85) 

P,^{k, E) = n^{k) 6{E- < E{k) >2Nc) (4.86) 
P,^{k, E) = niik) 5{E- < E{k) >3vc) • (4.87) 
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Figure 4.14: Different definitions of the scaling function vs. the scaling vari- 
able y. Dotted line: the longitudinal momentum distribution f l3.27p : solid 
line: the scaling function (14. 73 p . calculated with the spectral function; dashed 
line: the new scaling function f l4.77p . calculated with the nucleon momentum 
distributions. 
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Chapter 5 

Results of calculations of the 
inclusive cross section ratios 



Introduction 

In what follows, the results of our calculations of the inclusive cross section 
ratios shown in §1.2.31 will be presented. 



5.1 Inclusive cross section ratios in PWIA 

Let us recall that the experimental inclusive cross section ratio fll.4ip of 
nucleus A to the nucleus ^He was defined in §1.2.31 as follows 

which, within our novel approach, reduces to the following ratio 

-(A^ He) = ■ (5-2) 

Zaep + Naen FJjiy q) 

Our preliminary results of the PWIA ratio, for A = 56 are shown in Fig. 
15. Ij it exhibits a good agreement with CLAS data only for 1.5 ^ xbj ^ 2, 
i.e. in the region of 2NC; on the contrary, at xbj ^ 1-5 the PWIA does not 
lead to satisfactory results. 

This fact agrees with the results already shown in Fig. 14.81 in the region 
of 2NC the data of heavy nuclei scale to the Deuteron ones, and thus FSI 
effects vanish in the ratio r{A/^He), leading to the first plateaux; indeed in 
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Figure 5.1: The experimental cross section ratio shown in Fig. II. 141 compared 
with our preliminary PWIA theoretical results, for A = 56. 

the 2NC region, Eq. (15. 2 p reads as follows 

r(A, He) ~ — — - — — ^tt^ o„ , ^ ~ = const . (5.3) 

^ ' Zaep + Naen QH^^kdknl^'^ik) C^"- ^ ' 

In the kinematical region at xbj ^ 1-5, on the contrary, the ratio exhibits 
a strong sensitivity upon the A dependent FSI of the knocked nucleon with 
the residual system, and this is the reason of the disagreement between the 
experimental ratio and our calculations performed in PWIA. 
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5.2 FSI and distorted nucleon momentum dis- 
tribution 

In order to include FSI effects within our new model for inclusive scattering, 
we will now introduce the distorted nucleon momentum distributions, defined 
as follows [93] 

e) = j dv dr' e*P'n-(r-r') p^(r, r') . (5.4) 

The quantity appearing in Eq. (15. 4p . is the missing momentum 

Pm = q-k (5.5) 

defined in terms of the three-momentum transfer q and the momentum of 
the knocked out nucleon k (see §1.2.21) . with 6 being the angle between them; 
P£)(r,r') is the distorted one-body mixed density matrix, i.e. the quantity 

[94] 

. M <0o|Ft S^p{r,T')SF\<Po> 
pD{r,r) = ^ , I , ^ • (5.6) 

< (pom > 

In Eq. (15. 6p . 0o is the mean field wave function, = F 0o is the real- 
istic correlated wave function, S is the operator which takes into account 
FSI, and, eventually, /)(r, r') is the one-body density matrix operator. When 
5 = 1, PWIA is recovered, and the missing momentum equals the nucleon 
momentum before interaction. 

The distorted momentum distributions have been calculated in Ref. [M] 
adopting the eikonal Glauber representation for the quantity S, namely 

A 

S = l[ G{hi-h„zi-z,) (5.7) 
i=2 

where 

^(bi - hj, zi - Zj))l - e{zj - zi) r(bi - hj) (5.8) 
with r(bi — hj) being the usual Glauber profile function, i.e. 



r(b, - b,) = ^ e-^V^g . (5.9) 



Here, r = {bi,2;i} is the coordinate of the struck nucleon with transverse 
and longitudinal coordinates bi and zi (the axis z is along the direction of 
the struck nucleon), a^f^ is the total nucleon-nucleon cross section, bo is the 
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slope parameter of the total NN elastic cross section, and 6{zj — zi) ensures 
that the struck nucleon interacts only in the forward direction. 
Calculations similar of the ones of Ref. [91] for ^ife have been performed 
for closed shell nuclei ^^O and ^''Ca [22] • The results for the three nuclei are 
presented in Figs. 15.21 and IHT^ where they are compared with the distorted 
momentum distributions of the Deuteron. It is amazing to see that even in 
the case of the distorted momentum distribution a sort of Deuteron scaling 
is observed at high values of the missing momentum, where 

nj,{p^,e)^C^nj^{p^,e). (5.10) 

5.3 Inclusive cross section ratios with FSI 

Let us replace in Eq. (14.771) the undistorted nucleon momentum distributions 
with the distorted ones; we obtain 



where 



foiy) = foD{y) + f2D{y2) + hD{y^) (5.ii) 



Ip™I^Ip™I j <D{\vJ.o)de (5.12) 



y\ 

+ 00 

A 



f,^^{y,) = 2n \pjd\pj nU\PrnlO)de (5.13) 

J\y2\ J 

f3Diy3)=^ \Pm\d\Pm\ / ^^D ( I Pm I > ^) • (5-14) 

J\ys\ J 

By including FSI effects, the preliminary results shown in Fig. 15.41 are ob- 
tained, showing a good agreement with CLAS data, both at XBj ^ 1-5 and 
in the 2NC region (1.5 ^ xbj ^2), where Eq. (15.31) reads now as 

r{A,^ He) ~ = const . (5.15) 



5.4 3NC nucleon momentum distributions 

In order to include 3NC effects, and to extend our calculations in the region 
of 2 < XBj ^ 3, we need the 3NC nucleon momentum distributions n^{k), 
which are, to date, completely unknown. So the problem arises of how to 
determine them. Herebelow the following suggestion is presented. 
In Fig. the exact neutron spectral function of ^He and the one 
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calculated within the FNC model (12.641) are compared. The original idea of 
our approach is to subtract from the exact spectral function the FNC one, 
obtaining, by this way, the component of the spectral function due to 3NC, 
i.e. P3Ncik,E). 

The corresponding nucleon momentum distributions are nothing but the in- 
tegral of the spectral function 

r+oo 

ni{k)=47T / P,j,c{k,E)dE (5.16) 
so that the green area shown in Fig. 15.61 yields n^{k). 

Calculations of 3NC effects in the region 2 < x^j < 3 are in still in progress, 
and will be reported elsewhere [72], but the results of our preliminary calcu- 
lations show that 

n^{k)^n^{k) (5.17) 

which is in agreements with the finding summarized in Fig J1.15[ where the 
independent low momentum shell model nucleons and the 2NC high momen- 
tum nucleons almost exhaust the description of nucleons in a nucleus A. 
Moreover, in order to explain the second plateaux appearing in the inclusive 
cross section ratios (I1.14p . we should expect that n^{k) shows a ^He scaling 
behavior, i.e. 

ni{k)^Cinl"\k) (5.18) 



so that 



in the 3NC region, at 2 < x^j < 3. 
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(a) 



e=90 
e=o° 
e=i80° 




3 rr -U 4 



(c) 



Figure 5.2: The distorted nucleon momentum distributions of Eq. (15 ■4p 
versus the momentum |Pm,|, for (a) ^H, (b) ^He and (c) ^^O. Dashed hne: 
e = 90°; dot-dashed red line: 9 = 0°; dotted line: 9 = 180°; solid line: no 
FSI taken into account. After Ref. |92|. 
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e = i8o 

4. 




U 4 



(a) 



0=180 

16. 




3 . 4 



(b) 



Figure 5.3: Dashed line: the distorted nucleon momentum distribution of 
^He and ^^O at 6 = 180°; dotted line: the rescaled Deuteron momentum 
distribution at the same angle; full line: the undistorted momentum distri- 
bution. The same behavior appears at = 0° and 6 = 90°. After Ref. 

M. 
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Figure 5.4: The same as in Fig. I5.H with theoretical calculations performed 
with Eq. f l5TT]l . After Ref. |17|. 
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E [MeV] 

Figure 5.5: The neutron spectral function for ^He vs. the removal energy E, 
for two different fixed values of the momentum k. 
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E' [Me VI 



Figure 5.6: The 3NC contribution to the spectral function, extracted by 
subtracting from the exact spectral function the FNC one, both shown in 
Fig. 15. 5[ plotted versus the energy E* = E - E^^;. 
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Summary and conclusions 



In the first part of this Thesis, we have illustrated the theoretical techniques 
used to solve the many-body problem of nuclei and, in particular, the theoret- 
ical problems encountered in finding a solution at high values of the nucleon 
momentum. Prom a theoretical point of view, we have justified an approx- 
imation based only upon two- and three-body potentials, with higher order 
potentials being neglected. We have then pointed out the necessity of taking 
into account short range nucleon nucleon correlations (SRC) for a complete 
description of nuclei, and the necessity of a deep knowledge of them for an- 
swering different questions in various fields of the modern physics. 
In the second and original part of this Thesis, we have focused on the way 
SRC manifest themselves in the high momentum components of one- and 
two-nucleon momentum distributions; moreover, we have stressed out that, 
due to SRC, and in particular to the dominance of n-p pairs in a nucleus 
A, the nucleon momentum distribution for a nucleus A is nothing but the 
Deuteron momentum distribution rescaled by a constant C"^, a behavior 
known as Deuteron Scaling. 

We have then illustrated our novel approach to the study of SRC effects 
in inclusive lepton scattering off nuclei, based upon the introduction of the 
proper scaling variables y, y2 and y^; these take effectively into account the 
excitation energy E\_^ of the residual system in different ways, allowing us 
to describe the A(e, e')X cross section in terms of the corresponding momen- 
tum distributions generated by 2NC and 3NC. 

We have shown that, in the region of 2NC, our new approach, in terms of 
nucleon momentum distributions, results to be a good approximation of the 
one based upon the use of the spectral function. Moreover, we have illus- 
trated, that the experimental scaling function, in the 2NC region, not only 
scales to the Deuteron scaling function, but also exhibits A independent final 
state interaction (FSI) effects, mostly due to the FSI in the correlated pair. 
Calculations in the region of 3NC requires the knowledge of the nucleon mo- 
mentum distributions nf{k), and we have proposed an approach to obtain 
the three-nucleon momentum distributions from the knowledge of the full 
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5.4. 3NC nucleon momentum distributions 



and the 2NC spectral function. 

To sum up, we have demonstrated that the usual approach to scaling, based 
on the variable is not very useful as far as the investigation of the short 
range structure of nuclei is concerned. On the contrary, the direct, global, 
and A-independent link between the scaling function 7/2) and the lon- 

gitudinal momentum distributions f^{y) allows one to obtain information on 
the general behavior of the high momentum components in nuclei, which are 
governed by SRC. 

Eventually, we have demonstrated that, within our novel approach, the in- 
clusive cross section ratio r{A/A') reduces to the scaling function ratio of 
nuclei A and A'. Our preliminary calculations of the scaling function ratio, 
performed for A = 56 and A' = 3, show in PWIA a good agreement with 
CLAS data only for 1.5 < xbj ^ 2, i.e. in the region of 2NC; on the contrary, 
at XBj ^ 1-5, the PWIA does not lead to satisfactory results. This is not 
a surprising result, indeed, in the region of 2NC the data of heavy nuclei 
scale to the Deuteron ones, and thus FSI effects vanish in the ratio r{A/A'), 
leading to the first plateaux; in the kinematical region at xbj ^ 1-5, on the 
opposite, the ratio exhibits a strong sensitivity upon the A-dependent FSI 
of the nucleon knocked out from mean field states with the residual system. 
Including explicitly these FSI effects in the mean field contribution, we ob- 
tained a good reproduction of the experimental plateaux attributed to 2NC. 
Calculations including the 3NC configurations, which are necessary in order 
to extend our comparison with the CLAS experimental data to the region 
2 < Xfij < 3, are in progress, and will be presented elsewhere [72] . 
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